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We consider the neutrino (and antineutrino) flavors arriving at Earth for neutrinos produced in 
the annihilation of weakly interacting massive particles (WIMPs) in the Sun's core. Solar-matter 
effects on the flavor propagation of the resulting > GeV neutrinos are studied analytically within a 
density-matrix formalism. Matter effects, including mass-state level-crossings, influence the flavor 
fluxes considerably. The exposition herein is somewhat pedagogical, in that it starts with adiabatic 
evolution of single flavors from the Sun's center, with #13 set to zero, and progresses to fully realistic 
processing of the flavor ratios expected in WIMP decay, from the Sun's core to the Earth. In the 
fully realistic calculation, non-adiabatic level-crossing is included, as are possible nonzero values 
for 813 and the CP-violating phase 8. Due to resonance enhancement in matter, nonzero values of 
013 even smaller than a degree can noticeably affect flavor propagation. Both normal and inverted 
neutrino-mass hierarchies are considered. Our main conclusion is that measuring flavor ratios (in 
addition to energy spectra) of > GeV solar neutrinos can provide discrinination between WIMP 
models. In particular, we demonstrate the flavor differences at Earth for neutrinos from the two 
main classes of WIMP final states, namely W + W~ and 95% bb + 5% t + t~ . Conversely, if WIMP 
properties were to be learned from production in future accelerators, then the flavor ratios of > GeV 
solar neutrinos might be useful for inferring #13 and the mass hierarchy. From the full calculations, 
we find (and prove) some general features: a flavor-democratic flux produced at the Sun's core 
arrives at earth still flavor-democratic; for maximal #32 but arbitrary 821 and 813, the replacement 
8 — > n — 8 leaves the e flavor spectra unaltered but interchanges [i and r spectra at Earth; and, only 
for neutrinos in the inverted hierarchy and antineutrinos in the normal hierarchy is the dependence 
on the mixing phase 8 not suppressed to order Sm^i/ '8m%2- 

PACS numbers: 14.60.Pq, 95.85.Ry, 26.65.+t, 95.35.+d 



I. INTRODUCTION 

One of the most fundamental questions about our Uni- 
verse is its matter content. Recent cosmological obser- 
vations of type la supernovae [1], the cosmic microwave 
background [2] , galaxy-cluster evolution [3] , and gravita- 
tional lensing [4] show that known particle species con- 
tribute only Qb — 4% to the total energy density of 
the Cosmos. Abundance measurements of light elements 
together with the theory of Big-Bang Nucleosynthesis 
largely support these findings [5] . The above results also 
indicate the presence of two — as yet unknown — matter 
components: "dark energy" and "dark matter" with rel- 
ative densities of SIde — 73% and Qdm — 23%, respec- 
tively. Unraveling the nature of these novel types of mat- 
ter lies currently on the forefront of fundamental-physics 
research. 

One of the new components, dark matter (DM), ex- 
hibits gravitational clustering and is therefore believed 
to be composed of novel massive particles. Further com- 
pelling support for this idea is provided by theoreti- 
cal approaches to physics beyond the Standard Model, 
which typically require additional particles for consis- 
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tency. Weakly interacting massive particles (WIMPs) are 
a general and popular candidate in this context. A partic- 
ularly attractive example of a WIMP is the lightest super- 
symmetric particle (LSP), usually the neutralino, present 
in various SUSY extensions of the Standard Model [6]. 
This particle is a mixture of the four superpartners of 
the two neutral Higgs particles and the two neutral elec- 
troweak gauge bosons. The neutralino is expected to 
have a mass of order 10 2 GeV. The attractiveness of iden- 
tifying the LSP with dark matter is twofold: First, the 
LSP is stable, assuming i?-parity is an unbroken sym- 
metry. i?-parity is introduced to stabilize the proton 
to the level required by experiment. Second, a stan- 
dard thermodynamic calculation of decoupling yields the 
LSP abundance to be C>(lO _34 cm 2 /0Vi), where a a is the 
WIMP annihilation cross section for Majorana particles 
like the neutralino (for Dirac WIMPS, there is an addi- 
tional factor of Tdecoupie/AfwiMP ~ 1/30). Inputting the 
usual weak cross section apSs = -^wimp) ~ ^fI^i one 
naturally explains the inferred value Odm — 23% [7, 8]. 

Besides direct detection via WIMP scattering in cryo- 
genic detectors, indirect searches for WIMP annihilation 
into gamma rays, antimatter, and neutrinos are currently 
being pursued as promising experimental avenues [9, 10]. 
The gravitational fields of the Sun and the Earth can 
capture large numbers of WIMPs, which are expected 
to then infall toward the center of the body [11]. Since 
the annihilation rate is proportional to rip M , annihila- 
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tions of WIMPs captured by the Sun or Earth may be 
considerable. An important factor in WIMP capture is 
the surface area of the capturing body. Another is the 
efficiency of capture e cptr , which is related to the depth 
of the potential well relative to the mean kinetic energy 
of the infalling WIMP. The capture rate is proportional 
to the square of the body area and to the capture effi- 
ciency. There results a relative Sun-to-Earth capture rate 
of (e c s p u t n r /4 a t r r th ) (^o/^e) 2 , where R Q = 6.96 x 10 10 cm 
is the radius of the Sun and i?© = 6.38 x 10 s cm the 
radius of the Earth. Allowing for the l/D 2 flux dilu- 
tion between the source and the earthly detector, one 
then arrives at a naive relative flux ratio at Earth of 
F Su „/F Earth ~ (e s c ™/e*£ h ) (i? /i? e ) 2 (i?e/AU) 2 = 2 x 
10~ 5 (e^p t n r /e^p t r r th ). Three further pieces of physics more 
than compensate the 10 -5 factor and make the Sun the 
better source for experimental study. First, the capture 
efficiency of the Sun far exceeds that of the Earth, with 
relative potentials <~ Mj R that are 3000 times larger for 
the Sun. The second is that the solar potential reduces 
the lifetime of WIMP orbits among the planets, thereby 
reducing the capture rate at Earth [7]. The third factor 
is that center-pointing gravitational force ~ M(< r)/r 2 
is far greater for the Sun, due to its roughly exponential 
relation between density and radius. In fact, the DM 
in the Earth is not concentrated in Earth's center but 
rather broadly distributed in r. As a consequence, rip M 
is not optimized in the Earth, and the annihilation signal 
is fairly diffuse because it comes from a large fraction of 
a steradian. On the other hand, the Sun has a high cen- 
tral DM density, and an annihilation signal coming from 
a fraction of a square degree of solid angle. Thus, we will 
not consider further the WIMP annihilation flux com- 
ing from the Earth, but rather focus on the annihilation 
signal from the Sun. 

The solar medium absorbs all annihilation products 
except neutrinos with energies below tens of GeV. The 
typical energies of fusion neutrinos are 10 McV or less, 
while the typical energies of neutrinos and antineutrinos 
from the WIMP decay chain are GeV to tens of GeV. It 
follows then that the detection of higher-energy neutri- 
nos from the Sun provides an excellent tool for indirect 
WIMP detection. 

The idea of inferring WIMPs by measuring the associ- 
ated high-energy neutrino flux from the Sun's center has 
been the focus of numerous theoretical studies [12], and 
first experimental results have placed loose constraints on 
such scenarios [13]. Note, however, that neutrino detec- 
tion is highly dependent on neutrino flavor (v e versus 
versus v T ). These flavors evolve as the neutrinos travel in 
space and time. This is particularly so in the Sun, where 
the solar environment adds significant matter-dependent 
effects to the flavor evolution [14]. In fact, matter ef- 
fects do not only alter the flavor evolution, but they also 
treat the neutrino and antineutrino differently. There- 
fore, dark-matter searches via annihilation neutrinos typ- 
ically require considerable knowledge of neutrino mixing 
and mass parameters. 



Matter effects were not included in the early works 
cited in [12]. More recently, the MSW effect has been 
incorporated into the analysis of solar WIMP annihila- 
tion neutrinos by a number of authors [15-18]. Of these, 
only [18] is recent enough to have focused on the "last- 
solar-model-standing," the large mixing angle (LMA) so- 
lution. The other papers hedged their bets among LMA, 
SMA, LOW, and VAC solutions appropriate for their 
time when less was known [19]. 

Although our knowledge of neutrino parameters has 
advanced quickly, it remains incomplete. The neutrino 
mixing angles seem to converge to a "tribimaximal" 
form [20] with values in vacuum given to a good approx- 
imation by sin 2 623 ~ 1/2, sin 2 612 ~ 1/3, and sin 2 #13 <~ 
0. The large value of the "solar" angle #12 lies at the 
heart of the LMA solution. The neutrino mass-squared 
differences have converged to Sm^ <~ 2.5 x 10~ 3 cV 2 and 
8m\ 2 ~ 0.8 x 10~ 4 cV 2 . 

In Ref. [18], numerically precise terrestrial spectra re- 
sulting from the primary WIMP-annihilation channels 
are predicted, signal topologies in neutrino telescopes are 
given, and the possibility of discriminating among dark- 
matter models is discussed. Our investigation comple- 
ments that in Ref. [18] mainly in two ways. 

First, we focus on particularly promising observables: 
there has been a growing awareness in recent years that 
neutrino telescopes can infer flavor ratios of arriving 
neutrinos. It has been shown [21] that neutrino fla- 
vor ratios may be determined in ice and water neutrino 
telescopes from the observation of the ratio of muon- 
track events to jet events, where the former arise from 
charged-current (CC) scattering, and the latter from 
v e and v T CC events and all neutral-current (NC) events. 
The muon mean-free-path (MFP) before decay is cr M <~ 
6.3km(£' jU /GeV), which allows even 100 MeV muons to 
be identified. Large-volume magnetized-iron calorime- 
ters, liquid- Argon detectors, and active-scintillator detec- 
tors have the potential to perform even better neutrino- 
flavor identification [22]. The r MFP is ct t ~ 0.49 x 
10~ 2 cm (E T /GeV), leading to possibly identifiable t 
tracks for energies above 10 GeV. 

Second, we use a related but different treatment of the 
propagation of E <; GeV neutrinos from the solar core 
to Earth. We employ a density-matrix formalism and 
arrive at compact and convenient analytic results which 
we believe improve our intuition. 

Besides presenting matter effects that alter the flavor 
content of a propagating neutrino flux, possibly reso- 
nantly, solar matter also absorbs neutrinos above some 
energy i? a b s , due to the neutrino inelastic CC and NC 
interactions. We wish to neglect absorption by looking 
at neutrinos with E < E^. So we must estimate -E a bs- 
In the Sun, the number density of electrons n e decreases 
nearly exponentially with the distance r from the Sun's 
center. An approximate expression for the solar density 
profile n e (r) is given by [23] 

n e (r) = 245 N A £ - r/x& cm" 3 , A = 0.095 R Q , (1) 
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for r < Rq; Na is Avogadro's number. Equation (1) 
also serves as a rough estimate of the proton and neutron 
densities. 

Integration of this expression from r = to r = R Q 
gives the column density Sj from the Sun's center to its 
surface: 

S e = S p ~ S„ ~ 1.6 N A x 10 12 cm" 2 . (2) 

(For R 3> A Q , one also has Sj ~ rij(0) X & .) The neu- 
trino's optical depth r„ is the number of scattering MFPs 
from the Sun's center to the surface (the most probable 
number of interactions in transit); the fraction of pro- 
duced neutrinos that escape the Sun is e~ T " . The optical 
depth is equal to the product of cross section o v and col- 
umn density. Since ov./v ^ 0V e and S p + S n ~ 2 S e , we 
have for the optical depth t v ~ 2 <j v n S e . If we ask that 
a fraction / or more of the neutrino flux escapes from the 
Sun, we require 

T v <\nf-\ (3) 

With the onset of deep inelastic scattering at ~ GeV en- 
ergies, the neutrino and antineutrino cross sections grow 
linearly with E until TF-propagator effects become im- 
portant near E ~ tum ~ TeV. To a good approx- 
imation, the expression for the inelastic CC+NC neu- 
trino and antineutrino cross sections at energies between 
1 GeV and 1 TeV is [24]: 

<r vN (E) ~ 2 a pN (E) ~ 10- 38 cm 2 . (4) 

The rules of thumb arc that the neutrino cross section 
is twice the antineutrino cross-section, and the CC con- 
tribution is a bit over twice the NC contribution. Also, 
the fractional energy loss per scatter is about 1/2 for 
neutrinos at the energies relevant here, and half that 
again for antineutrinos. In CC scattering, the neutrinos 
and antineutrinos are lost upon the first scatter; in NC 
scattering, the neutrinos and antineutrinos persevere but 
with the fractional energy losses just presented. Since the 
NC/CC scattering ratio is less than a half, we may and 
will neglect NC scattering in what follows. 

The energy E a ^ s can now be calculated by substituting 
Eq. (4) and twice times Eq. (2) for r„ in Eq. (3). The 
result is -E a bs = 50 (In/ -1 ) GeV for neutrinos, and twice 
that for antineutrinos. Some explicit numerical results 

for Kbs = 2 Kbs are 10 GeV when / = 90% : 36 GeV 
when / = 70%, and 70 GeV when / = 50%. If we 

discount the NC processing of neutrinos and attribute 

absorption to only the CC, then E a b s is 50% larger. A 

careful numerical calculation of absorption versus energy 

is given in Fig. 1 of Cirelli et al. [18]. The conclusion there 

and here is that up to energies of ~ 10 GeV, absorption in 

the Sun can be neglected, and above 100 GeV absorption 

is severe. 

The format of this paper is as follows. Flavor oscilla- 
tions are presented in Sec. II. Matter effects, both adi- 
abatic and non-adiabatic, the two possible mass hierar- 
chies, and zero versus nonzero #13 and 8 are discussed in 



detail. In particular, two interesting theorems for the 5 
dependence of the neutrino flavors, valid for the adiabatic 
regime, are presented. The application of the formalism 
to the solar WIMP annihilation is made in Sec. III. A 
summary and conclusions comprise Sec. IV. Some tech- 
nical details are collected in four appendices. 



II. PHYSICS OF FLAVOR MIXING 

This section presents a density-matrix approach for 
the flavor evolution of high-energy neutrinos injected at 
the center of the Sun. In the first two subsections, cer- 
tain simplifying assumptions regarding neutrino-mixing 
parameters and neutrino propagation in matter are em- 
ployed; our goal is to set up our formalism and notation, 
to provide intuition about the involved effects, and to 
set the stage for a more sophisticated study in the third 
subsection. 



A. Review: oscillations without matter 

In the absence of a matter background (i.e., in vac- 
uum), the free- neutrino Hamiltonian is diagonal in the 
mass basis: 

H M =£l + -^diag(m 2 ,m 2 ,m 2 )+0(m 4 /£ 3 ). (5) 

Here, E is the neutrino's energy, and 1 is the 3x3 unit 
matrix. We label the mass eigenvectors in vacuum as 
\j) = |1), 1 2), 1 3). However, neutrinos are produced by 
the weak interaction. In this production basis, called the 
"weak basis" or "flavor basis," the interaction eigenstates 
are labeled by the flavors of their charged-current part- 
ners, as I a) = |e), \fi), |r). In general, the flavor basis 
is rotated with respect to the mass basis, which leads to 
periodic flavor oscillations among propagating neutrinos. 
The oscillation length is given by 

4tt£ 

v = Mk 

-•Ho*)^)*- <6 > 

where Sm 2 k = m 2 — m| denotes the neutrino mass- 
squared differences. As mentioned in Sec. I, the values 
of 5m 2 - k inferred from solar and atmospheric neutrino- 
oscillation experiments are Sm^i ~ 10~ 4 cV 2 and 5m§ 2 ~ 
10~ 3 eV 2 , respectively. Thus, the vacuum oscillation 
lengths are well contained within the Sun for energies 
below 30 GeV. 

The transformation between the mass basis and the 
flavor basis is described by the vacuum mixing matrix 
U a j = (j\oi) = with a = e, jj,, t and j — 1, 2, 3. 

The conventional parametrization of the vacuum mixing 
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matrix invokes three angles and a phase [25]: 

U = i? 23 (M Uj R 13 (6 13 ) Us R 12 (8i 2 ) , (7) 

where Rjk(8jk) determines a rotation in the jk plane by 
an angle 9jk- The phase matrix 



U S = 



e lS / 2 
1 

e- lS / 2 



(8) 



may 



be thought of as complexifying R\ 3 {9i 3 ) to 



R 13 (6 13 ,6) = ulR 13 Us 

cos 9i 3 



er l5 sin 6»i3 

1 



-e +iS sin 6>i3 cos#i 3 



(9) 



We have omitted two additional Majorana phases, as 
they leave unaffected neutrino oscillations. 

As mentioned before, the central values of the mix- 
ing angles inferred from oscillation experiments are quite 
consistent with the tribimaximal angles [20] given explic- 
itly by 6*32 = 45°, 9 2 i = 35.26 . . .°, and 6 13 = 0. The full 
tribimaximal mixing matrix is 



U = 



( 



\ 







(10) 



We remind the reader that rows are labeled from top to 
bottom by the flavor indices e, /x, and r, and columns 
are labeled left to right by mass-eigenstate indices 1, 2, 
and 3. 

Due to decoherence of the propagating neutrino 
phases, we shall see that it is the classical probabilities 
\U a j | 2 , obtained by squaring the mixing-matrix elements, 
that are of interest to us. We collect these squared ele- 
ments into a matrix that we denote as U. For tribimax- 
imal mixing, its explicit form is 



(11) 



One cannot avoid noticing the simplicity of the tribimax- 
imal probabilities: the v e contents of v\, u 2 , and v 3 are 
2/3, 1/3, and 0, respectively, with the remaining proba- 
bility being shared equally between and v T . 

The 3-<7 experimental limits on the mixing-angle 
ranges are [26] 




0.34 < sin 2 9 32 < 0.68 , 



0.24 < sin 2 6*2i < 0.41 , 



sin 2 6>i3 < 0.044, or 8 13 < 12° 



(12) 



(13) 



(14) 



Notice that with U e3 = sin 6*13, the "zero" in U may 
be as large as 4%. Later in this paper we will explore 
the sensitivity of our results to experimentally allowed 
variations in 9\ 3 . 

Neglecting matter effects, the density matrix for neu- 
trino production in the Sun is 



w e |e)(e| + Wf, +w T \t)(t\ 



(15) 



where the weights w a are the relative flavor ratios at 
production normalized to obey w e + w M + w T = 1. The 
density matrix for antincutrino production is given by the 
same formula when the w a are replaced with antincutrino 
production weights w a ■ The first and second expressions 
in Eq. (15) present the production density matrix in the 
flavor and mass bases, respectively. The point of exhibit- 
ing p in the mass basis is to twofold. First, this basis 
diagonalizes the free-neutrino Hamiltonian and is there- 
fore appropriate for describing neutrino propagation to 
Earth. Second, it facilitates a discussion of decoherence, 
as is explained next. 

As the neutrinos propagate from the Sun to the Earth, 
each off-diagonal operator \j)(k\ acquires a phase factor 
exp(— i<fikj) w hh large phase <f>kj = L8m 2 .^/2E 1 where L 
is the distance from neutrino-production site to the de- 
tection site at Earth. The size of the Sun's DM core and 
the detector's size contribute a large 5L, and energy res- 
olution in the detector contributes a SE to the variance 
in the phase. For a statistical sample of events, these un- 
certainties randomize the phases, effectively equating the 
phase factor to its mean value ( cxp(— i<j)ki) ) = 0. Thus, 
off-diagonal elements of p, associated with the initial co- 
herence of the neutrino wave function, may be dropped. 
The loss of phase information allows us to work with 
classical probabilities. We effectively obtain 



p(Earth) = J2 w a \U aj \ 2 \j)(j\ 



(16) 



for the density matrix describing neutrinos arriving at 
Earth. 

With the above result, the flavor probabilities at Earth 
are given by 



(/3|p(Earth)|/3) = ]>> Q \U aj \ 2 \U 0J \ 2 . (17) 



In matrix form, this equation becomes 













Z e 1 









(18) 



Inputting the tribimaximal values leads to the explicit 
"flavor-propagation matrix" 



UU 1 = 




(19) 
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and to 

P VQ ^ Ve = [lOiUe + 4(«v +w T )} 

= l - (2 + 3<), (20) 

= 1(7-3^). (21) 

The equations for antineutrino flavor probabilities are ob- 
tained from these neutrino results via the replacement 
w a — > Wa- 
in Eq. (21), the two probabilities P„ _>„ and P Vq - >Vt 
are equal. Also, in the first expressions on the right 
in Eqs. (20) and (21), and w T enter symmetrically. 
These results are manifestations of a v^-Vr interchange 
symmetry inherent in the tribimaximal mixing matrix. 
The v^-v T interchange symmetry is exact for #23 = 45° 
(maximal v il -v T mixing) and #13 = [27]. The inter- 
change symmetry holds in matter as well as in vacuum 
because the v^-v T sector remains unaffected by the mat- 
ter potential. More manifestations of the interchange 
symmetry will become evident in expressions and fig- 
ures to be presented in subsequent sections. In the 
second expressions on the right in Eqs. (20) and (21), 
w e + + w T = 1 has been implemented. It is because 
of the v^-Vt interchange symmetry that the final result 
depends on the single weight w e . This greatly reduces 
the latitude in characterizing the flavor distribution at 
the source. 

For neutrinos with energies well below any matter reso- 
nances, the matter effects are negligible. These neutrinos 
therefore propagate as if in vacuum, and for tribimaxi- 
mal mixing they produce flavor probabilities at Earth 
given by the formulas above. At higher energies, matter 
effects become substantial. One test for the efficacy of 
the solar matter potential on neutrino mixing would be 
to determine flavor probabilities with values other than 
those presented in Eqs. (20) and (21). 



B. Oscillations with matter: 
adiabatic approximation 

In matter, the forward scattering amplitude of the neu- 
trino leads to an effective potential [14]. The scatter- 
ing amplitudes common to all neutrino flavors generate 
a common potential, which is unobservable in neutrino 
oscillations. However, there is a scattering amplitude 
unique to v e and v e flavors. It is the VF-exchange diagram 
for scattering on background electrons, with exchanges in 
the t-channel for v e and in the s-channel for v e . Thus, 
there results an effective potential for the v e and v e fla- 
vor states, given by V e (r) — \f2Gp n e (r) for v e , and — V e 
for u e , where n e is the electron number density. The ef- 
fect of the matter contribution on neutrino flavor mixing 



is best found by first transforming the vacuum Hamilto- 
nian (5) to the flavor basis, labeled by {|e), \ |r)}, and 
then adding to this the flavor-diagonal effective potential. 
Then, the flavor-basis effective Hamiltonian in matter 

Hp = u diag(-*m| 0, 6ml 2 ) ^ + ^ ^ 

emerges. In Eq. (22), some contributions to Hp propor- 
tional to the identity have been omitted for convenience 
because they leave unchanged the flavor mixing, as ex- 
plained earlier. For antincutrinos, the sign of matter- 
potential term is reversed and U is replaced by U* . 

The matrix that diagonalizes the flavor Hamiltonian in 
the matter background is called U m : Hm = Hp U m . 
Such a matrix is only fixed modulo the ordering and the 
phases of its column vectors. The phases are without 
physical significance because they can be absorbed into 
the definition of the states | j, r) defined below. The or- 
dering of the column vectors, on the other hand, deter- 
mines the order of the eigenvalues along the diagonal of 
Hm- Our convention is this: the order of the column 
vectors in U m is such that the diagonal entries of Hm 
exhibit the same ordering in magnitude as those of the 
matrix M — diag(— Sm^i, 0, 5m\ 2 ) in Eq. (22). Note that 
this convention depends on the mass hierarchy, via the 
sign of 6m^ 2 - This convention is chosen such that U m 
can smoothly approach the vacuum mixing matrix U in 
the limit of vanishing matter potential. 

Given a position-dependent matter potential, the mix- 
ing matrix U m (r) will also depend on position. For the 
matrix which diagonalizes Hp at the Sun 's center, we re- 
serve the notation U m without argument r. In the solar 
core, the matter density plateaus, and the matter poten- 
tial at the Sun's center is 

V e (0) = 70V2Gf AUckT 3 = 5.6 x 10~ 12 cV. (23) 

Thus, neutrinos from the solar core with energy above ~ 
Sm 2 jk /2V e (0) - ((5m2 fe /10- 5 cV 2 )MeV will feel' the effect 
of the Sun's matter potential. This class includes the 
E <; GeV neutrinos from WIMP annihilation. 

With the mixing matrix in matter U m (r) and the fla- 
vor eigenstates at hand, one can construct the "instan- 
taneous" eigenstates of Hp via \j,r) = J7^(r)|a). Al- 
though these states are strictly speaking not solutions 
to the equations of motion, they will turn out to be in- 
valuable in what follows. For future reference, we fur- 
ther define the set of "matter states" {|l) m , |2) m , |3) m } 
as the exact solution of the Schrodinger Equation with 
the Hamiltonian Hp given in Eq. (22). It is typically dif- 
ficult to determine explicit expressions for these states. 
In this section, we simply note that equating | j) m to \ j, r) 
is known as the "adiabatic approximation." We defer to 
Sec. II B 4 a more detailed discussion of this approxima- 
tion. To summarize, the weak-interaction Hamiltonian 
is diagonal in the flavor basis |a). On the other hand, 
the free-neutrino Hamiltonian in vacuum or matter is di- 
agonal in the mass basis \j) or the matter basis \j) m , 
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respectively. An approximation for the matter states is 
given by the instantaneous cigcnstates \j, r) . 



1. Diagonalization of Hf for 9n = 

The diagonalization of Hp is simplified by the obser- 
vational input that #13 is small, so we can initially ap- 
proximate it by zero. Then, the 3x3 vacuum mixing 
matrix simplifies to U = ^23(^23) ^12(^12)- This permits 
us to cast Hp into block-diagonal form via the transfor- 
mation H' F = i?23 1 (^23) Hf #23(023) because the rotation 
matrix i? 2 3 commutes with the matter-potential term in 
Eq. (22). The result is 



H. 



2x2 










2E 



(24) 



with 



H 



2x2 



V e (r) 



Sm 21 2 a 

^COS e 12 



5m. 



21 



sin 2#i2 



V 



-^f sm20 12 



5m%i o 

— sin 2 1 

2E 



'12 



(25) 



determining the upper 2x2 block. Note that the matter 
potential V e (r) and off-diagonal elements are confined to 
H 2x2 . This means that the |3) m mass eigenstate in mat- 
ter is unaffected by background electrons and decouples 
from the other two states. This allows a two-state treat- 
ment of the original three-flavor problem leading to some 
simplifications, which we discuss next. 

Since the state |3) m is unaffected by the background 
matter, it is identical to the vacuum mass eigenstate 
I 3 ) = {KxhaWa) = sin6»23|Ai) +cos6» 23 |t). The fact 
that |3) m = 1 3) decouples from the other two states im- 
plies that the currently unknown neutrino mass hierarchy 
(777,3 > m 2 > mi vs. 77i2 > mi > 7773) is unimportant in 
the present case characterized by #13 = 0. 

Another simplification following from the three- to 
two-state reduction is the elimination of any possible CP- 
violation from the effective mixing matrix. This is evi- 
dent in that the phase 5 does not enter into Eqs. (24) and 
(25). Still, the neutrino and antineutrino cases have to 
be treated separately, for their respective Hamiltonians 
differ in the sign of the matter-potential term. 

A third simplifying feature arising from the reduction 
to two-state mixing is the straightforward recognition of a 
resonance condition. Mixing becomes resonant (maximal 
angle and degenerate eigenvalues) when the two diagonal 
elements oi H 2x2 become equal, i.e., when 



E R {r) =Sm 2 21 cos 29 12 /2 V e {r) 



(26) 



The resonant energy at the Sun's center is therefore 

E R (0) = cos 26»i2 (SmlJl CT 4 cV 2 ) f MeV . (27) 



The tribimaximal value for 008 26*12 is 1/3. Equation (1) 
implies that outside the solar core the matter poten- 
tial decreases, so that the resonance energy increases. 
Thus, only neutrinos with energies exceeding the solar- 
core value (27) can experience resonant conversion. The 
sign of Sm 21 is fixed as positive by solar data. It follows 
that for antineutrinos the right-hand side of Eq. (26) is 
negative, and there is no possibility of resonance conver- 
sion when 6*i3 = 0. 

The next step is to determine the eigenstates |l) m and 
|2) m of the effective 2x2 Hamiltonian (25). Follow- 
ing the conventional analysis [29], we approximate these 
states by the instantaneous eigenvectors \l,r) and |2,r) 
that diagonalize #2x2 at the point r. This adiabatic ap- 
proximation is valid if the matter density changes slowly 
enough with distance. In this section, we assume adia- 
baticity holds. In the next section, we relax this assump- 
tion. 

In the limit 



<5m 2 /2V;(0) 



(28) 



the instantaneous eigenstate |1,0) at the Sun's core r ~ 
is predominantly determined by the matter potential, 
leading to |1,0) ~ |e). The majority of neutrinos from 
solar WIMP annihilation are expected to have energies 
E <; O(GcV); for these neutrinos, |1,0) ~ |e) is a good 
approximation [28]. Then the only mixing is due to 6*23 
and we have U m = ^3(623) [28]. Thus, the remaining 
instantaneous mass state is given by |2, 0) ~ cos 6*23 — 
sin 6*23 1 t). 

As for the vacuum case in the previous subsection, we 
employ a density-matrix analysis to predict the flavor 
evolution. We again begin with the density matrix p(0) 
describing the mixed-flavor ensemble produced at the so- 
lar core. Our above considerations regarding the instan- 
taneous mass states at the Sun's center yield [28] 



p(0) = w e |e)(e| +Wp |m> <M| + w t \t)(t\ 

-|3) + |2,0)\ /(3| + (2,0| 



= w e |l,0)<l,0|+u;„ 
+ w 



V2 J\ V2 
|3)-|2,0)\/<3|-(2,0r 



V2 



V2 



(29) 



In the latter expression, the best-fit value of #23 — 45° 
has been implemented. 

We continue by studying the propagation of the neu- 
trinos from the Sun's core to its surface. Quantum in- 
terference effects can be safely neglected because the 
neutrinos quickly decohere after production: With the 



exception of possible resonance points 



where the 



two eigenvalues of H 2x2 approach each other, the os- 
cillation lengths A are small compared to solar scales. 
For example, in the production region r ~ 0, the 
set of approximate instantaneous eigenvalues of H' F 
is {-<5T77 2 2 /2£;,(5777 2 2 /2i;, V e (0)}. The relative sizes of 
(57773 2 1 IE and V e (0) dominate, so that the oscillation 
scale is given by E/8m\ 2 ^ O(100 km) for E > 1 GeV, 
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5m 2 2 w 10~ 3 eV 2 , and V^ -1 (0) — 40 km. If the uncer- 
tainty in SL is comparable to these oscillation lengths, 
then the phase of the oscillation is randomized per neu- 
trino and averaged out in the statistical sense. Such is 
the case here, since the size ~ 0.01i? Q -\/100 GeV/mDM 
of the WIMP annihilation region in the Sun's core [11] 
exceeds the oscillation scale of 100 km. Due to this de- 
coherence, we can transform the density matrix p to the 
mass basis, and then disregard the off-diagonal pieces, as 
was done for the vacuum case. The result is three in- 
coherent fluxes associated with the matter states 1 1 , 0) , 
1 2, 0), and |3) emerging from the production region at the 
center of the Sun. 

Neutrinos corresponding to |3) will remain in this state 
because it is an exact eigenstate of the Hamiltonian un- 
affected by the solar electron density. However, the in- 
stantaneous eigenstates 1 1, r) and |2,r) change with re- 
spect to the flavor basis as r increases, and transitions 
between these states can occur. The adiabatic approxi- 
mation neglects these transitions, so that a particle will 
remain in an instantaneous eigenstate. We note that the 
flavor-oscillation wavelengths for GeV neutrinos close to 
resonances in the Sun are long compared to the 0.1 Rq 
scale length of the exponentially decreasing solar density, 
so that the adiabatic approximation breaks down in some 
circumstances. The adiabatic analysis presented here is 
nevertheless interesting, for it is instructive, simple, and 
applicable in certain energy ranges, as we shall see. 

We are now in the position to establish the evolution 
of the density matrix. The discussion above shows that 
at sufficiently large neutrino energies the vacuum masses 
are negligible compared to the effective mass of the state 
1 1 , 0) . However, outside the Sun, where V e (r > Rq) = 0, 
the mass hierarchy is reversed: the instantaneous states 
are just the vacuum mass states 1 1, r > Rq) = |1) and 
|2,r > Rq) = |2) characterized by mi < m 2 . Since the 
usual quantum-mechanical level repulsion prohibits level 
crossing, the states in the adiabatic approximation evolve 
as follows: as the neutrino travels to the solar surface, we 
have |1,0) -> |2> and |2,0) -> |1) up to phases [28]. As a 
result, we obtain the density matrix 

p(r > Rq) = w e |2)(2| + i (w, + Wt ) (|3)(3| + |1)(1|) (30) 

for neutrinos emerging from the Sun. Upon exiting 
the Sun, the \j) states become exact eigenvectors of 
the neutrino Hamiltonian, so that transitions between 
these states cease to occur. Outside the Sun in empty 
space, the density matrix therefore remains effectively 
unchanged obeying p(Earth) = p(r > Rq). 

This result can be employed to determine the relative 
flavor fluxes measured in Earth-based experiments. The 
probability of detecting the neutrino flavor (3 <G {e,u.,r} 
by such experiments is 

= (/3|p(Earth)|/3) 

- w e \U p2 \ 2 + \ K + w r) (1 - \U ( n\ 2 ) ■ (31) 



In the final step, we have used Sj=i W[ij\ 2 — 1, which 
arises from the unitarity of U. The v^-v-t interchange 
symmetry is evident again. The relation We+w^+Wr = 1 
allows us to cast the expression in terms of w e only. We 
remind the reader that although this result incorporates 
solar matter effects, the \Upj\ 2 in the final expression are 
elements of the vacuum matrix U. 

Note that an amusing result emerges if \U e2 \ 2 = 
|t^»2| 2 = | 1 2 , namely, that all three flavors arrive at 
Earth with equal probability 1/3, independent of the fla- 
vor ratios at production. While \U e 2\ 2 = \U^ 2 \ 2 = \U T 2\ 2 
may seem contrived, in fact it is the cornerstone of the 
tribiximal mixing matrix. We discuss this next. 



2. Neutrino flavor evolution 

To obtain the expected neutrino flavor probabilities at 
Earth, we insert the tribimaximal mixing values for U 
given in Eq. (11), into Eq. (31). Doing so, and using the 
normalization relation ^2 a w a = 1, we arrive at 

Pvq— *V e = PfQ = Pvq— *V T = g • (32) 

As advertised above, this result for neutrinos is indepen- 
dent of the w a which characterize the flavor distribution 
at the Sun's center. Such a result could hardly be sim- 
pler (or more democratic). It is easily understood as 
well: With the assumption that #13 = 0, the state |3) 
is completely decoupled, and so the higher-energy reso- 
nance is absent. Thus, a v F produced with energy above 
the lower-energy resonance emerges from the Sun as pure 
\v 2 ) in the assumed adiabatic limit. With tribimaximal 
mixing, the decomposition of \v 2 ) into flavors gives the 
observed ratios \U e2 \ 2 : p^] 2 : \U t2 \ 2 = 1 : 1 : 1. On the 
other hand, z/^'s and i/ T 's produced in the Sun's center 
equilibrate with each other, a consequence of the v lx -v T 
interchange symmetry. They then emerge from the Sun 
with equal parts of the two states orthogonal to v 2 , i.e., 
with equal parts of v\ and ^3 . The flavor decomposition 
of these states is |t//3i| 2 + IC//33I 2 = 1 - |£^/32| 2 , which is 
again 1:1:1 for the tribimaximal mixing matrix. Thus, 
v e and independently and v T generate a 1 : 1 : 1 flavor 
ratio at Earth. As a consequence, the 1:1:1 predic- 
tion is independent of the initial w a . The essence of this 
argument is visible in Eqs. (30) and (31). 

Of course, the state evolution is not perfectly adia- 
batic, and 613 is probably not exactly zero. With either 
alteration, one expects deviation from the 1:1:1 fla- 
vor ratios. We address these issues in Sec. II C. Fig- 
ure 2 shows a comparison of our simplified approach de- 
scribed here with the full non-adiabatic result obtained 
in Sec. IIC. 

We remark that the democracy derived above for the 
arriving neutrino flavors is a nontrivial prediction. First 
of all, it depends on the approximate validity of the tribi- 
maximal mixing matrix. Second, it depends on the solar 
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matter, as outlined above. Reference to the vacuum pre- 
dictions given in Eqs. (20) and (21) shows that without 
the solar-matter effect, flavor democracy results only for 
the single value u e — 1/3. Third, even with tribimaxi- 
mal mixing and the solar-matter effect, democracy is not 
a prediction for the arriving antineutrino flavors, as we 
show next. 



3. Antineutrino flavor evolution 

Antincutrinos from the Sun yield more interesting fla- 
vor probabilities. The construction of these flavor proba- 
bilities follows the same logic as those for neutrinos. How- 
ever, there is one major difference: the hierarchy of the 
effective mass in the solar core and the vacuum mass hi- 
erarchy are identical, so that no mass-label permutation 
is necessary. To obtain the antineutrino analogue of Eq. 
(30), we must therefore take |1,0) -> |1) and |2,0) -> |2) 
in Eq. (29). This gives 

- ^e|^l| 2 + ^K+^)(l-|t//3l| 2 ). (33) 

The tribimaximal mixing probabilities are \U e \ | 2 = | and 
\U^\ 2 = \U Tl \ 2 = \. It follows that 

P Po _ Pe = ^ (1 + 3 w g ), 

P^, = Pu^f> T = ^ (5-3 w e ) . (34) 

Here, the answer docs depend on the flavor ratios w & at 
production in the Sun's center. For example, an injection 
flux of purely D e leads to P Po ^ Pc = | and Pv^^v^ = 
Py Q ^y T = whereas an injection flux with no v e leads 
to Pp Q ^p e = \ and -Pp ^ Pfl = Pp ->i? T = j2- Thus, in 
principle the w & are measurable. Their determination 
would help unravel the nature of the DM source. 

Why are the antineutrino flavor ratios not democratic, 
as were the neutrino ratios? The answer lies in the fact 
that v e emerges as pure v\ , which is distributed in flavor 
as 4 : 1 : 1, and the P M and D T emerge as equal parts 
v-i and P3, which are collectively distributed as 2 : 5 : 5. 
Thus, the resulting flavor ratios depend on the initial 
ratios w & . 

Since resonances in the antineutrino sector are ab- 
sent within the normal hierarchy, another question arises. 
Why do the antineutrino results here differ from the vac- 
uum results of Section II A? The difference is due to the 
mixing effect at the antineutrino production site. The 
large matter potential causes |e) to be almost purely |l) m , 
whereas in vacuum |e) would be distributed among the 
mass states \j) as |f/ e j| 2 - 

In Fig. 2, the approximate antineutrino probabilities 
derived here for adiabatic evolution and #13 = can be 
compared to the full analysis performed in the next sec- 
tion. It is seen in Fig. 2 that the simple approxima- 
tion characterized by constant (i.e., energy- independent) 



probabilities, is good over a relevant range of several 
decades in energy for the antineutrino. 

We end this section with the remark that in both the 
neutrino and the antineutrino cases, all flavor probabil- 
ities converge to 1/3 (democratic 1:1:1) for w e and 
w s = 1/3. This is because for these latter values all mass 
eigenstates are equally represented at production, and 
subsequent level repulsion simply permutes the labels of 
these mass states. This result holds for any values of 
the mixing matrix, and it also applies in the presence of 
matter effects, as we prove in Sec. II C 2. 



4. Validity range of the above simplified results 

As is evident in Fig. 2, for the normal mass hierarchy, 
the adiabatic approximation is good for neutrinos with 
energies between ~ 10 MeV and ~ 10 GeV, and for an- 
tineutrinos with energies above ~ 10 MeV. What physics 
sets these limits of validity? 

The low-energy limit of validity E\ ow is determined by 
the requirement that at production, the matter state |e) 
be well separated from the other states; in turn, this 
requires that the matter-potential term dominates the 
vacuum mass term in the effective Hamiltonian (22) at 
production. A sufficient condition for this stipulation is 
given by Eq. (28), which yields E\ ow ~ 70 MeV when the 
value of 5m\ 2 is taken. In fact, the limit E\ ow can even be 
lower due to the particular form of the mixing matrices 
in the Hamiltonian (22). 

At energies well below the lowest resonance, the mat- 
ter term in the Hamiltonian becomes negligible and vac- 
uum results are approached. As can be seen in Fig. 2, 
the low-energy asymptotes for all neutrino and antineu- 
trino curves have the vacuum values given in Eqs. (20) 
and (21). 

The high-energy limit of validity -Ehigh is determined 
by the onset of non-adiabaticity, i.e., when appreciable 
transitions between the instantaneous eigenstates begin 
to occur. The non-adiabaticity scale Ena is established 
in the next section. Borrowing results from the next sec- 
tion, we can conclude the following: Ena — > 00 for the 
nonresonant v case, and from Eq. (51), Ena ~ 100 GeV 
for the resonant v case with the parameters listed in 
Fig. 2. 

Neutrino absorption in the Sun is significant at energies 
above ^100 GeV, thus providing another limit on the 
validity of our analysis. If the absorption is universal 
across the flavors, then one might expect a reduction in 
flux and a change in spectrum, but little or no alteration 
of the flavor results we derive in this paper. 



5. Neutrinos and antineutrinos together 

An unmagnctized detector is ill-suited to resolve v 
events from v events. For such a detector it is more 
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appropriate to consider the weighted average 

0*1/ Pvq P^Q — >vp ~t~ @i> Pi 



(P. 



P- 



o v F v& + crp Fi 



(35) 



as an observable. Here, F„ B and F Pq are the neutrino 
and antineutrino fluxes produced by WIMP annihilation 
in the Sun, and o v and ct p are the neutrino and antineu- 



J 



trino cross sections in the detector. We define the neu- 
trino to antineutrino cross-section ratio as r a = Up/a,,. 
Experimentally, r a ~ 1/2 in the energy range of inter- 
est. In situations in which neutrino production is charge 
symmetric, the initial v a and v a fluxes are equal, as are 
w a and Wo- The weighted average of Eqs. (31) and (33) 
is then simply given by 



(P 



VQ^vp at up 



) = 



[^e (3 1*7/32 | 2 



l) + l-\U P 2\ 2 ] + r a [w e (3\U 01 \ 2 



l) + l-l^ll 2 ] 



2(1 + r a ) 



(36) 



r 



Assuming the tribimaximal mixing angles, one obtains 

4 + r CT (5-3w e ) 



(-Pi/0 —>!/,, or 
(Pisq — >v e or v, 



v^l — {Pv Q ^v T orv T ) — 

2 + r ff (3w e + l) 



12(l + r CT ) 



6(l + rv) 



(37) 



The w e dependence of these results, although weakened 
compared to the pure v sample given in Eq. (34) , implies 
that detectors that sum v and v events still have some 
discriminatory power to resolve different source models. 
In Sec. Ill we discuss this possibility. 



6. Summary of adiabatic evolution in matter 

In our description of adiabatic evolution in matter, 
we adopted three simplifying assumptions governing neu- 
trino (and antineutrino) production by solar WIMP an- 
nihilation, and subsequent propagation out of the Sun 
to Earth. These assumptions are that all neutrino en- 
ergies exceed the resonant energy at the Sun's center 
(-Er(O) ~ 10 MeV), that the neutrinos evolve adiabat- 
ically through the Sun, and that neutrino mixing angles 
are given by the tribimaximal values, including 613 = 0. 
As a consequence, we obtained energy-independent fla- 
vor probabilities for the neutrinos arriving at Earth. We 
learned that solar matter plays a significant role in deter- 
mining the expected flavor ratios at Earth. However, the 
higher-energy resonant value has yet to play a role, 
since #13 = decouples the |3) state from the |1) and 
1 2) states. In the next subsection, we relax the various 
assumptions. 



C. Non-adiabatic oscillations and nonzero O13 



1. Nonzero 613 and three-flavor mixing 

We begin again with the density matrix p„(0) at pro- 
duction, expressed in flavor space: 



/°"(°) = ^2 W <* \ a )( a \ 



(38) 



where a <E {e,fi,r} runs over the three flavors, and the 
w a are the relative flavor fluxes at production, as before. 
Paralleling our previous simplified analysis, we expand 
the flavor cigenstates \a) in the basis of instantaneous 
eigenstates \ j, r) at r = 0, where j € {1, 2, 3} labels again 
the effective-mass eigenvalue: \a) — \j, 0). 

We obtain 



MO) = w ° 



(u m ) aj (u m y ak |j,o)(fc,o| 



(39) 



Recall the earlier definition of U m as the solar-core in- 
stantaneous mixing matrix in matter, i.e., (U m ) a k = 
(a\k,0). The matter at the Sun's center characterizes 
the core environment where the neutrinos are produced. 
Again, as so defined, U m diagonalizcs Hp only at the 
center of the Sun. 

The same decoherence arguments presented in the pre- 
vious section allow us again to neglect the off-diagonal 
density-matrix elements in Eq. (39). Therefore, in what 
follows we will take 



^(0)-^«; Q |([/ m ) aj | 2 |j,0)0-,0| 



(40) 



as the neutrino-production density-matrix. As a check, 
we note that with a large matter potential isolating the 
\v e ) state and maximal #32 = 45°, Eq. (40) becomes iden- 
tical to our earlier Eq. (29). 



The physics becomes more complicated when simpli- 
fying assumptions are abandoned. For example, if #13 
is no longer assumed to be zero, then the state |3) no 
longer decouples and a full three-flavor analysis with two 
resonant energies is typically required. 



2. Non-adiabatic evolution 

Next, we re-examine the evolution of this system dur- 
ing propagation to the solar surface, but this time with- 
out the adiabatic assumption of the previous section. 
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That is, we now allow for transitions between the in- 
stantaneous eigenstates \j, r) . 

Once the neutrinos reach the vacuum at the Sun's sur- 
face, transitions between the eigenstates \j,r > i? Q ) = 
\j) no longer occur because in free space the \ j) are exact 
eigenstates of the Hamiltonian. For our analysis, we are 
therefore interested in the solar-propagation amplitudes 
Akj = (k\j, 0) . Embedded in these amplitudes are the im- 
portant non-adiabatic transitions between instantaneous 
eigenstates, most likely to occur near resonance. 

In terms of the amplitudes A k j, the density matrix 
given in Eq. (40) evolves to become, beyond the Sun's 
surface, 

Pv {r >R Q )= E w « \{Um) aj \ 2 A* u A kj \k){l\ . (41) 

a,j,k,l 

It is seen here that non-adiabatic transitions replenish 
the off-diagonal elements of p. However, the A k j in gen- 
eral contain position-dependent phases, which must be 
averaged over the production region. Furthermore, ad- 
ditional loss of phase information occurs as the neutri- 
nos continue en route to Earth, as discussed in Sec. II A. 
So again the off-diagonal terms average to zero. Con- 
sequently, Eq. (41) is reduced to the more manageable 
expression 

p„(Earth) = ]T w a \(U m ) aj \ 2 \A kj \ 2 \k)(k\ . (42) 

The amplitude Aj k is seen to enter the final formula only 
via the jump probability P k j = \A k j -\ 2 , as expected in the 
absence of coherence. 

Equation (42) can now be employed to determine the 
probability of detecting a neutrino of flavor j3 at Earth. 
We obtain 

(/3|p(Earth)|/3) 
= E w a \(U m ) aj \ 2 P kj \U 0k \ 2 . (43) 

a >3>k 

In matrix form, this equation is given by 

/ Pvq^u. \ ( We \ 

P*e^ \=UPUl\wA. (44) 

V Pv @ ^Vr J \ W r J 

Here, we have defined the matter matrix (U m ) a j = 
\(U m ) a j\ , in analogy to our earlier vacuum definition 
Equations (43) and (44) generalize Eqs. (31) and (18) in 
the previous subsections. In the adiabatic (P k j — > 5 k j) 
and matter-free (U m — > U) limits, these equations reduce 
to their analogues in the previous subsections. 

With Eq. (43) at hand, we may establish rigorously 
an interesting feature of the flavor evolution: If the Sa- 
vors at production are democratically distributed (w e — 
= w T = 1/3 ), then the flavor probabilities observed 
at Earth will also be democratically distributed. In other 
words, the particular ratio (3,3,3) of flavors is left un- 
changed by the combination of all oscillation and matter 



effects. This is expected in that off-diagonal density- 
matrix elements average to zero, so that the density- 
matrix evolves by changing bases, from flavor in matter 
to evolving mass in matter, to constant mass in vacuum, 
and finally to flavor at earth; but a density matrix pro- 
portional to the unit matrix initially, remains so through- 
out basis changes. 

We verify this conclusion by setting w a = 1/3 in 
Eq. (43), and rewriting the remaining terms slightly: 

p »^»> = l E i^i 2 E p ^ E Wm) aj \ 2 ■ (45) 

k j a 

The unitarity of U m implies that the sum over a appear- 
ing on the far right-hand side of Eq. (45) equals unity. 
Similarly, the sum of elements in the row or column of 
P k j must add up to 1 because of probability conserva- 
tion. The remaining sum over k is also equal to 1 due 
to the unitarity of U. We are thus left with the claimed 
result that P VQ ^ vp = 3 for all flavors (3. Consequently, 
any experimentally inferred flavor ratios deviating from 
1:1:1 exclude a democratic source. 

Equation (43) allows predictions to be made for neu- 
trino flavors at Earth, given an initial distribution of fla- 
vors in the solar core. It has been used extensively in the 
context of suppression of the expected solar v e rate at 
Earth. For neutrinos from the solar fusion cycle, one has 
w e = 1, with all other w a and w & equal to zero. How- 
ever, little study of Eq. (43) has been performed for an 
ensemble of solar flavors, as would arise from WIMP an- 
nihilation in the Sun. An early study of neutrino flavor 
physics arising from solar WIMP annihilation was per- 
formed in Ref. [15]. Much has been learned about the 
neutrino mixing matrix in the intervening fifteen years. 
We incorporate subsequent knowledge into our work here. 

An approximation for U m is derived in Appendix A, 
so our remaining task is to present explicit expressions 
for the transition probabilities P k j = \A k j\ 2 . In the adia- 
batic approximation of the previous section, these |^4fej | 2 
were just 6 k j. However, the adiabatic approximation can 
break down in resonance regions, where the oscillation 
length increases by the factor 1 / sin 29 in the two-state 
approximation. If the oscillation length becomes compa- 
rable to or larger than the scale length of the matter po- 
tential I din V e (r)/dr| _1 , the medium and therefore the 
matter-induced mixing angles are changing sufficiently 
fast, so that adiabaticity is violated [30]. When neutrino 
evolution is non-adiabatic, P k j has a complicated struc- 
ture and includes off-diagonal pieces. The nature of P k j 
is discussed next. 

For further progress, we employ another piece of ex- 
perimental data: dm 2 ,^ and 5m 2 2 differ by more than 
an order of magnitude, so that both resonances in our 
three-neutrino system can be treated separately and in- 
volve only two levels at a time [31]. For each of the two 
resonances, an improved Landau-Zener approximation 
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establishes a variant of the Parke formula [30] 

P c = e(E-E Rm exp(-rsin^) -cxp(-r) 

1 - exp(— r) 

as the generic leading-order form of the level-crossing 
probability [29]. Here, 9 is the effective two-level vac- 
uum mixing angle, Y is the "adiabaticity parameter", 
and Er(0) = Er is the resonance energy at the so- 
lar core. These quantities depend upon the three-level 
vacuum neutrino parameters, the matter-density profile, 
and the resonance under consideration. Before applying 
Eq. (46) to each resonance in the present situation, we 
discuss some of its general features. 

As discussed in Sec. II B 1, outside the Sun's center the 
resonance energy increases Er(t ^ 0) > Er- It follows 
that only neutrinos produced with energies E > Er can 
experience resonance conversion on their way to the Sun's 
surface. Since a resonance is required for appreciable 
level-crossing probabilities, we need E > Er in order 
for P c 7^ 0. This feature is ensured by the presence of 
the threshold function 9 in Eq. (46). We also note that 
the crossing probability P c goes smoothly to 100% as 
9 approaches zero. This is expected and necessary; it 
reflects the fact that two states cannot repel each other 
when they are decoupled. 

Explicit expressions for the adiabaticity parameter T 
show [29] that it depends on the neutrino energy as l/E. 
Thus, it is useful to define a "non-adiabatic energy" Ena 
via the simple equation 



r = 



Ena 
~E~ 



(47) 



Ena sets the energy scale for non-adiabatic effects. In 
the limit E — > 0, Y — > oo and the crossing probability 
P c goes to zero. This is the adiabatic limit. The onset 
of non-adiabaticity is found by expanding Eq. (46) under 
the condition E <C Ena- One has 



P c = cxp(-rsin 2 9) - exp(-r) + 



(48) 



Taking P c = e~ 3 ~5%asa characteristic value for the 
onset of non-adiabaticity, the leading term in Eq. (48) 
determines an onset energy Ena for non-adiabatic effects: 



Ena = 2 sin 2 9 E NA ■ 



(49) 



Thus we have an energy condition for nonadiabaticity: 
E > max{i?R, Ena}- Note, however, that for very large 
energy E >• Ena, Y goes to zero and P c approaches 
cos 2 9. This reflects the fact that at very high energies 
the oscillation wavelength exceeds the matter region, and 
so vacuum oscillations result. 

For large sin 2 9, the definition (49) may be inconsistent 
because the next-to-leading contribution in Eq. (48) can 
become sizable. However, we will see below that only 
the cases 9 = 9\i , 6*13 < 45° are of interest in the present 
context. The ratio of the second to first terms in the 



expansion (48) evaluated at Ena is exp(— 3 cot 2 9). Even 
for the ultra-conservative assumption 9 — 45°, this ratio 
is less than 5%, so that Eq. (49) is more than sufficient 
for our purposes. 

Although the generic crossing probability (46) provides 
the basis for the determination of Pjk, many details also 
depend on the number of resonances encountered (0, 1, 
or 2), the mass hierarchy, and on whether neutrinos or 
antineutrinos are considered. For example, with #13 = 
the state |3) decouples as we have seen, and there is but a 
single resonance; it lies in the neutrino sector with either 
the normal or inverted mass hierarchies. With #13 7^ 0, 
there are two resonances; they both lie in the neutrino 
sector with the normal mass hierarchy, but lie one each in 
the neutrino and antineutrino sectors with the inverted 
mass hierarchy. The schematic Fig. 1 displays the possi- 
bilities. We use superscripts I and h on relevant variables 
to denote the lower- and higher-energy resonances, and 
+ and — to denote the normal and inverted hierarchies, 
respectively. Some obvious energy regions where evolu- 
tion is purely adiabatic are (i) E v < E l NA , (ii) all E„ 
with the normal hierarchy, (iii) all Ep with the inverted 
hierarchy when #13 = 0, and (iv) E D < E NA with the 
inverted hierarchy when #13 7^ 0. We will consider below 
the various cases, normal vs. inverted hierarchy, and #13 
zero vs. nonzero. 

In principle, a neutrino produced on the back side of 
the Sun but traveling toward us may encounter two res- 
onances before leaving the production region. There is a 
standard formalism available to describe this physics [29] . 
However, we may ignore this "res-in-res-out" possibility 
for two reasons. The first is that the region of DM an- 
nihilation is expected to be so small (^ 0.01 Rq) that 
only neutrinos within a particular narrow energy band 
could encounter the resonances. The second reason is 
that the matter density within 0.1 Rq of the solar core 
deviates from the exponential behavior in Eq. (1): it be- 
comes nearly constant [23] . This near-constancy of the 
density profile holds even better within Adm ~ 0.01 Rq. 



3. (Anti) Neutrino flavors - normal hierarchy 

For the normal hierarchy 5m\ 2 > 0, there are two res- 
onances for neutrinos and none for antineutrinos. The 
lower-energy neutrino resonance occurs between the two 
lightest states. For this situation, the characteristic 
quantities Er, Ena = YE, and 9 in the crossing- 
probability formula (46) are given by [29, 32] 
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These expressions are valid at leading order. Appropriate 
to the case which we are discussing, the superscripts I 
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FIG. 1: Schematic illustration of the relative positions of the lower- and higher-energy resonances for the normal and inverted 
mass hierarchies. Neutrino (antineutrino) eigenvalues are on the right (left) of the vertical effective-mass axis. The dotted lines 
indicate non-adiabatic transitions. 



and + label the lower-energy resonance and the normal 
hierarchy with 8m\ 2 > 0> 

The logarithmic derivative inside the absolute-value 
sign is to be evaluated at the radial position tr of the 
respective resonance. In the present case of an exponen- 
tial density profile, this factor is equal to A© ~ 0.1 Rq 
and independent of r, which simplifies the analysis. In 
particular, we may set 



E l + A = 7T A 5m 2 21 ~ 110 GeV - "_ a 2 \ t2 ■ (51) 



10" 4 eV 



The corresponding expression for the onset of non- 



adiabaticity (49) in the present case is 
1 



E l+ 



■ sin 2 9 12 E l + A . 



(52) 



For the higher-energy neutrino resonance occurring be- 
tween the two heaviest states, the characteristic quanti- 
ties En, Ena, and 9 are given at leading order by [29, 33] 
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FIG. 2: Solar neutrino and antineutrino flavor probabilities at Earth versus energy, for a single injection flavor and for the 
normal mass hierarchy. The v e , v^, and v T spectra at Earth are shown as solid, gray, and dashed lines. The neutrino mass- 
squared differences are <5m|i = 8.0 x f(T 5 eV 2 and Smh = 3.0 x 10~ 3 eV 2 , and the mixing angles are #12 = 35° , 6*13 = 0° , and 
#23 = 45° . The CP- violating phase S is set to zero. The vertical dotted lines mark the characteristic scales for the lower-energy 
resonance given in Eqs. (50) and (52). There is no higher-energy resonance when #13 = 0. Note the v^-Vr and v^-Vt interchange 
symmetries discussed in the text. Because #13 is set to zero here, these results remain valid also for the inverted-hierarchy 
case. Note that the high-energy neutrino fluxes are identical to the low-energy antineutrino fluxes. This accidental feature is 
attributable to tribimaximal mixing, as explained in Appendix B. 



where the superscript h+ labels the present higher- 
energy resonance, normal-hierarchy situation. From 
Eq. (49), we obtain here 

E h N + = ±sm 2 e 13 E h N + A (54) 

for the onset of non-adiabaticity. 

We now assemble all the pieces necessary for includ- 
ing non-adiabatic matter effects into our analysis for the 
normal-hierarchy case. The neutrino level-crossing prob- 



abilities are determined by 

(\-P l c P[ 0\/l \ 

Phi - Pi l-P l c l-P* P c h , (55) 

V 1 y v p c !-Pc J 

where the superscripts I and h imply evaluation of the 
crossing probability (46) at the lower resonance with 
Eqs. (50) and at the higher resonance with Eqs. (53), re- 
spectively. We remark that with two resonances at play, 
Pjk Pkj hi general, so that the ordering of the two non- 
commuting matrices on the right-hand side of Eq. (55) is 
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FIG. 3: Solar neutrino and antineutrino flavor probabilities at Earth versus energy, for a single injection flavor and for the 
normal mass hierarchy. Here, we have taken #13 = 12°, 5 = 0. All other neutrino parameters are as in Fig. 2. The and v T 
spectra and and v T spectra are interchanged if 5 = tt is chosen. Vertical dotted lines mark the characteristic scales for the 
lower-energy resonance given by Eqs. (50) and (52) and the higher-energy resonance given by Eqs. (53) and (54). 



essential. 

On the other hand, for antineutrinos there are no res- 
onances present in the normal-hierarchy case. The adia- 
batic relation Pjy = 5k j is exact. 

The terrestrial flavor fluxes resulting from pure flavor 
injection in the Sun's center are plotted in Figs. 2 and 
3 for representative sets of neutrino parameters. One 
sees that the processed flavors arriving at Earth are far 
from pure. The neutrino spectra in Fig. 2 illustrates the 
three energy regions that may result from the influence 
of a matter resonance. At E v > E^ A , the flavor evo- 
lution is complicated by non-adiabaticity. In the region 
Br ^ E v < Ej^ A , the flavor ratios have the simple adia- 
batic values presented and discussed in Sec. II B 2. Below 



Eft, the matter state |e) is not well separated from the 
other states, and the description is again complicated. A 
simplifying feature for all energies is the \x-r degeneracy 
that necessarily results when #32 = 45° and #13 = 0. 

Figs. 2 and 3 differ in that the minimal #13 = is input 
into the first, and the 3-cr maximal #13 = 12° is input into 
the second. Several qualitative differences result when 
#13 is nonzero. The most noticeable is that with #13 ^ 0, 
the mixing between |3) and an injected |e) significantly 
suppresses the v e flavor at Earth in the adiabatic high- 
energy region between E R + and E^. The next most 
noticeable change is the breaking of the \x-t degeneracy, 
for both neutrinos and antineutrinos. As explained in 
Sec. IIC 6, the and v T spectra, and the and v T 
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spectra, are interchanged if 6 = is replaced with 5 = n. 

Although there are no antineutrino resonances with 
the normal hierarchy, there are nevertheless features in 
the antineutrino flavor spectrum at the resonant ener- 
gies. These "kincmatical reflections" arc due to the 
following physics: As the energy passes from below to 
above the resonant value, the effective Hamiltonian (22) 
morphs from being vacuum dominated to being matter 
dominated. In turn, the mixing angle(s) associated with 
the resonance morph from their vacuum value(s) to their 
matter value (s). If the respective signs of the vacuum and 
matter parts of the Hamiltonian arc opposite, as happens 
for neutrinos in the normal hierarchy, then the transition 
may be dramatic; at resonance the mixing angle passes 
through 45°. If instead, the signs of the two Hamiltonian 
pieces are the same, as with antincutrinos in the nor- 
mal hierarchy, then the transition in mixing-angles is less 
dramatic, but possibly still significant. The transition of 
mixing angles through the lower resonant energy in the 
v e and v T curves of Fig. 3 are evident. The top antineu- 
trino panel with pure v e injected is the simplest to under- 
stand. Below resonant energies, the vacuum angles are 
relevant, and we expect a flavor ratio of (5/9, 2/9, 2/9), 
as derived from Eq. (19) or cquivalcntly from Eqs. (20) 
and (21). Above the resonant energies, the injected \v e ) 
is identified in matter with the state |1). In vacuum 
(e.g., at Earth) and to lowest order in sin #13 cos S, this 
mass state decomposes into the flavor ratio (2/3, (1 + 
3V2 sin6»i 3 cos<5)/6, (1 - 3 \/2 sin6>i 3 cos5)/6). At the 
3-er maximal value of 3-\/2 sin #13 cos<5 = 0.90, the final 
v r fraction is driven almost to zero. 



of the 1-2 resonance involves neither |3, r) nor the value 
of 6m^2- The resulting jump probability for neutrinos is 



kj 



l-P l c P l c 

1 



(56) 



where P c is evaluated with Eqs. (46), (50), and the re- 
placement ?H ► I—. Note that |l,r) and |2,r) are now 

the two heaviest states. 

For antincutrinos, the resonance involves 1 1 , r) and 
1 3, r), which are the two lightest states in the inverted 
mass hierarchy. An argument paralleling that leading to 
Eq. (53) yields for this case 
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Here, the superscript h— denotes the higher-energy 
inverted-hierarchy resonance. The energy for the onset 
of non-adiabaticity is given by 



E NA - 3 sin 013 E NA ■ 



(58) 



and the jump-probability matrix now takes the form 



1-P C ' 1 P c h 
P kj =[ 1 

Pi 1 I -Pi 1 



(59) 



4- (Anti) Neutrino flavors - inverted hierarchy 

Next, we consider the inverted hierarchy, which is char- 
acterized by 8m\ 2 < 0- Because the (decoherent) fluxes 
are unaffected by the mass hierarchy when #13 = 0, as dis- 
cussed in Sec. II B 1, Fig. 2 remains valid for either mass 
hierarchy. In the more general situation with #13 7^ 0, it 
is no longer true that the |3) eigenstate of the Hamilto- 
nian completely decouples, and so a few additional con- 
siderations are necessary. First, note that only one of 
the two resonances occurs for neutrinos; the other one 
now appears in the antineutrino sector. In particular, 
the condition of a large resonance separation, which af- 
fects the quality of our approximation, is met trivially for 
the inverted hierarchy. Second, for the analysis of each 
resonance we can employ the same bases and approxi- 
mations as we used for the normal hierarchy because the 
analysis essentially depends on the experimental result 
that |<5m§ 2 | » l^m^l independent of the sign of <5m§ 2 . 

Armed with these considerations, we find that the reso- 
nance between the states |1, r) and \2,r) stays in the neu- 
trino sector, and the present results are equal to those in 
the lower-resonance normal-hierarchy situation. Hence, 

we only need to replace the superscripts l-\ ► I— in Eqs. 

(50) and (52). This result is expected since the analysis 



Here, P£ is determined by Eqs. (46) and (57). The 
inverted-hierarchy flavor fluxes on Earth for a represen- 
tative set of neutrino parameters are shown in Fig. 4. 
A comparison of Figs. 3 and 4 reveals that the effect of 
moving the higher-energy resonance from the neutrino to 
the antineutrino sector via the change in mass hierarchy 
is dramatic. 

The density-matrix formalism presented above accom- 
modates any mixture of flavor fluxes produced in the so- 
lar core. For illustration, we have selected only examples 
with pure single-flavor initial states in Figs. 2, 3, and 4. 
The previously discussed sources are effectively incoher- 
ent, which implies that the fluxes for mixed- flavor cases 
can be inferred from these figures by simply taking the 
appropriately weighted incoherent sum of our pure-flavor 
results. 



5. Sensitivity to O13 

Because #13 is experimentally bounded to be less than 
about 12°, its measurement will be more difficult than 
that of the other two mixing angles. It is therefore in- 
teresting to discuss the #13 sensitivity of the terrestrial 
flavor ratios of neutrinos from solar WIMP annihilation. 
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FIG. 4: Neutrino and antineutrino flavor probabilities on Earth versus energy, for the inverted hierarchy. Here, we have taken 
<5w-32 = —3.0 x 10 _3 eV 2 . All other neutrino parameters are as in Fig. 3 (including 613 = 12° and 5 = 0). The and v T 
spectra and and v T spectra are interchanged if 8 = 7v is chosen. 



Comparison of Fig. 2 with Figs. 3 and 4 reveals that 
the relative flavor fluxes on Earth indeed depend on 
813. The normal-hierarchy antineutrinos in Fig. 3 and 
the inverted-hierarchy neutrinos in Fig. 4 show changes 
with respect to Fig. 2 that are of the expected size 
sin 12° ~ 0.2. However, the changes for normal-hierarchy 
neutrinos and inverted-hierarchy antineutrinos are more 
significant. The question arises as to why a relatively 
small change in 613 can lead to effects of this magnitude. 

The answer to this question lies in the transition from 
adiabaticity to non-adiabaticity for the higher resonance, 
which is governed by #13: Equations (54) or (58) show 
that the onset of non-adiabaticity directly scales with 
sin 2 #13. In Fig. 2, we have taken #13 = 0. The loss of 
adiabaticity therefore occurs at zero energy. Moreover, 



the level-crossing probability (46) becomes Pj} = 1 for 
zero mixing angle. In other words, the |3) m state is de- 
coupled and does not participate in the adiabatic label 
permutation resulting from level repulsion, as discussed 
in Sec. II Bl. In Figs. 3 and 4, however, #13 7^ and 
|3) m is not decoupled. If in such a situation an adiabatic 
region exists — as in the case for sin #13 = 12° — the state 
experiences level permutation |3) m — > |2) m while prop- 
agating to the solar surface. It is this mechanism that 
causes the large 6*13 effects for normal-hierarchy neutri- 
nos and inverted-hierarchy antineutrinos. 

How small can #13 be for this effect to occur? Since 
the label permutation is valid only for energies above 
resonance in the adiabatic regime, we must require that 
E lies above the higher resonance and below the onset of 
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non-adiabaticity for this resonance, i.e., 

E h R +< E <E h NA , 

E h R < E < E h NA , (60) 

for normal-hierarchy neutrinos and inverted-hierarchy 
antineutrinos, respectively. Here, the various energy 
scales are given by Eqs. (53), (54), (57), and (58). It 
follows that the onset of non-adiabaticity must lie above 
the resonance energy. Such will be the case for values of 
#13 down to a fraction of a degree. 

We conclude that the terrestrial flavor fluxes from neu- 
trinos originating in solar WIMP decay are quite sensi- 
tive to the mixing angle #i 3 . Accordingly, the study of 
solar- WIMP annihilation properties would greatly bene- 
fit from a pre-knowledge of the size of #13. On the other 
hand, if dark-matter properties can be extracted from 
other experiments, such as direct WIMP discovery in un- 
derground detectors, or even better, at the LHC or ILC, 
then a measurement of #13 through WIMP annihilation 
in the Sun may be feasible. 

6. The CP-violating phase 5 

Currently, the CP-violating phase 8 is unconstrained 
by experiments. We are thus led to study the effects of 
8 on the propagation of high-energy solar neutrinos. 

The general non-adiabatic density-matrix framework 
discussed in Sec. II C 2 has made no assumptions regard- 
ing the size of the phase 5. The energy range of interest 
for solar neutrinos from WIMP decay is a few 100 MeV to 
about 100 GeV. Below this range, matter effects are sup- 
pressed and the standard vacuum results are approached. 
Above this range, absorption effects become dominant, 
and our framework is no longer accurate. Since this en- 
ergy range is mainly below E^ and E l ^ A , the adiabatic 
approximation is useful. 

We again take 6*12 and #32 close to their observed tribi- 
maximal values given in Sec. II A. For the remaining 
mixing angle, the largest experimentally allowed value 
#13 = 12° is most interesting: From Eq. (9) one may in- 
fer that in a CP-conserving calculation, as here for deco- 
hered neutrino propagation, the 6 dependence will always 
occur in the combination 3?(e l5 sin #13) = cos (5 sin #13. It 
is therefore only in combination with larger values of #13 
that one may seek observable 8 effects. Thus, when illus- 
trating effects of nonzero 8, we will generally dramatize 
the result by taking #13 to be 12°, its maximally allowed 
3-cr value. 

Concerning the <5-parameter, two important analytic 
features emerge. The first is Feature (i): For # 32 = 45°, 
but arbitrary 621 and 8i 3 , the replacement 8 — > n— 8 inter- 
changes the role of and v T , and and v T , and leaves 
unaffected the terrestrial v e and v e fluxes. Although some 
versions of this result are known, we have included a proof 
in Appendix C for completeness. We remark that Fea- 
ture (i) makes no assumptions about the neutrino energy 



or adiabaticity, so that its validity actually extends be- 
yond our energy range of interest, from a few 100 MeV 
to about 100 GeV. 

Assuming continuity, Feature (i) implies that for 8 = 
90°, the and v r fluxes are identical. This is consistent 
with the previously determined condition 5i(t/ e 3) = for 
v^-Vt interchange symmetry: inspection of the mixing 
matrix U shows that U e z = —i sin #13 is purely imaginary 
for 8 = 90°. 

The second result is Feature (ii): For # 32 and #12 
at their tribimaximal values, the terrestrial flavor ratios 
for normal-hierarchy neutrinos and inverted-hierarchy 
antineutrinos are independent of 8 to zeroth order in 
(Sm^i/Sm^)- In addition, the terrestrial v e and v e fluxes 
are independent of 5 to zeroth order in (Sm^i/Sm^) for 
any hierarchy. This feature holds only in a finite energy 
region determined by adiabaticity and by V e E » 8m 1 - k . 
Fortunately, this energy region coincides with our range 
of interest, from a few 100 MeV to about 100 GeV; and 
fortunately again, (8m\ x l 8m\ 2 ) ~ 1/30 is a small param- 
eter. We verify this result in Appendix D. 

As mentioned above, #13 = eliminates 8 from the 
mixing matrices U and U m . It follows that Fig. 2, for 
instance, is unaffected by the value of 8. Figures 3 and 4, 
on the other hand, should exhibit some changes when 8 is 
varied. With Features (i) and (ii) at hand, certain phase 
effects on the terrestrial flavor fluxes can be understood 
qualitatively. An immediate consequence of Features (i) 
and (ii) is the degeneracy of the and v T flavor spec- 
tra in the adiabatic region for normal-hierarchy neutrinos 
and inverted-hierarchy antineutrinos. This characteristic 
is apparent in Figs. 3 and 4. Feature (ii) implies that 8 
effects are confined to the v ^ and v T fluxes in Fig. 3 and 
to the and v T fluxes in Fig. 4. We must look at these 
fluxes when seeking any sizable 8 dependence. 

The qualitative behavior of these fluxes under changes 
of 8 is dictated by Feature (i). As 8 — > 90°, 8 and tt—8 ap- 
proach each other, the [i- and r-flavor fluxes must move 
toward each other; at <5 = 90° they become degener- 
ate. When 8 is increased further, the /x- and r-flavor 
fluxes reverse their role with respect to each other. The 
first-order results for A in Appendix D establish that 
the 8 dependence of these fluxes is of the generic form 
B + C sin #13 cos 8, where B and C depend on the mass 
hierarchy, on v versus P, and on #13. The precise expres- 
sions for B and C in each case can be inferred from Eqs. 
(D9) and (D13). 

In Fig. 5, the 8 dependences of the various fluxes for 
inverted hierarchy (IH) neutrinos and normal hierarchy 
(NH) antineutrinos are shown. The NH neutrino and 
IH antineutrino spectra are not shown since their cos 8 
dependence is suppressed by the factor (Sm^ / '8m\ 2 ) ~ 
1/30. Although the neutrino energy in the figure is taken 
to be E = 10 GeV, the results remain approximately un- 
changed as long as E stays within our selected energy 
range. Apparent in Fig. 5 arc the features discussed 
above. In particular, the cos 8 <-> — cos 8 equivalence to 
the [i <-> t flavor interchange is evident. 
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FIG. 5: Terrestrial neutrino fluxes versus the phase 8. The neutrino mass-squared differences are ^mJi = 8.0 x 1CT 5 eV 2 
and 5m%2 = ±3.0 x 10~ 3 eV 2 , where the upper (lower) sign refers to the normal (inverted) hierarchy. The mixing angles are 
#12 = 35°, #i3 = 12°, and #23 = 45°, and the neutrino energy is E = lOGeV. The /i- and r-flavor fluxes exhibit opposite cos 8 
dependences, whereas the i^ e -flavor fluxes are approximately constant. Curves for normal-hierarchy neutrinos and inverted- 
hierarchy antineutrinos are not shown, as they are all nearly constant as a function of 8. These various results illustrate Features 
(i) and (ii) derived in the text. 



III. WIMP-MODEL PREDICTIONS FOR 
EARTHLY NEUTRINO FLAVORS 

In general, the WIMP LSP is a superposition of four 
Majorana neutralino states. These are the two neutral 
Higgsinos (SUSY partners of the two Higgs doublets), the 
Bino (SUSY partner of the 5(7(2)-singlet vector bo- 
son), and the neutral Wino (SUSY partner of the neutral 
W° member of the 5'C/(2)-triplet gauge boson). WMAP 
data yield Hdm — 23%, which in turn constrains the de- 
coupling of the WIMP and therefore bounds the WIMP 
mass and superposition. The non-observation of WIMPs 
in Tcvatron data at Fermilab further constrains combi- 
nations of WIMP mass and superposition. However, the 
fact remains that almost any neutralino superposition 



can be made consistent with present limits [34]. Theo- 
retical arguments are largely expressions of preference or 
prejudice. For example, proponents of "gauge-mediated 
SUSY breaking" are led to a "Bino-like" LSP, while pro- 
ponents of "anomaly-mediated SUSY breaking" are led 
instead to a "Wino-like" LSP. Thus, there is little guid- 
ance as to what neutralino superposition comprises the 
LSP. 

Fortunately, two disjoint sets of annihilation channels 
are common in LSP models. One (I) is the W + W~ chan- 
nel. The other (II) is the bb + t + t~ channel. There are 
further channels available to WIMP annihilation beyond 
(I) and (II), but these are calculated to be subdominant. 
Higgsino and Wino annihilation favor channel I, if the 
LSP mass exceeds M\y; Higgsino annihilation but not 
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Wino annihilation also produces Z Z if above threshold, 
but the neutrino spectrum resulting from this channel 
is similar to that from W + W~ , so we need not con- 
sider Z Z separately. Bino annihilation favors channel 
II, although even a small amount of mixing of the Bino 
with other neutralinos is likely to lead to a dominance 
of W + W~ production. For an LSP with mass less than 
Mw, the decay channel of the LSP is (II). 

Since the annihilating WIMPS are a CP-symmetric ini- 
tial state, the neutrino and antineutrino spectra emerg- 
ing from WIMP annihilation must be identical. The fla- 
vor processing by solar matter, described in this work, is 
not CP-symmetric. Therefore, differences between neu- 
trino and antineutrino flavor ratios at Earth are expected. 
Conversely, any observed differences between neutrino 
and antineutrino flavor ratios at Earth are an indication 
of the solar matter effect. In principle, it is possible to 
separate neutrino and antineutrino spectra in proposed, 
large, magnetized iron-calorimeter experiments. 

It is not hard to understand the reasons for the two 
preferred decay channels. Pairs of S'?7(2)-doublet Higgsi- 
nos or Winos have an electroweak coupling to W + W~ 
suppressed only by the t-channel W exchange. How- 
ever, the 5{/(2)-singlet Bino does not have this cou- 
pling. On the other hand, when the W + W~ channel 
is not available, all LSP pairs, being Majorana parti- 
cles, annihilate to massless final-state pairs in an I = 1 
(p-wave) angular- momentum state suppressed by /? DM , 
where /3dm is the velocity of the LSP. With massive 
final-state pairs, chirality flips enable an s-wave decay 
suppressed by the factor /3dm ( m Lai/^DM)- The veloc- 
ity of the non-relativistic LSPs is sufficiently small that 
the s-wave dominates and one gets the 6 6 final state 
Ncml/(N c ml + m^) = 95% of the time, and the t+t~ 
final state m^./(Ncml + m*) = 5% of the time; here, 
Nc = 3 is the number of QCD colors. The neutrinos 
from tau decay r — > v T + • • • are more copious and pro- 
vide a harder spectrum than those from b decay, so it 
turns out that the 5% t+ t~ mode makes a significant 
contribution to the final neutrino flux. 

In the figures to follow, in which all w a are taken from 
Ref. [18], we will display the flavor ratios to be expected 
at Earth for the two WIMP annihilation classes (I) and 
(II). For further insight, we will also show separately 
results for the pure b b and pure t + t~ decay modes. 
However, before displaying the flavor ratios expected at 
Earth, we show in Fig. 6 the flavor ratios produced by 
WIMP annihilation in the Sun for each of the four de- 
cay modes. These flavor spectra are qualitatively easy 
to understand. The 6 6 production and subsequent decay 
to (charm+) t+v t is phase-space suppressed relative to 
decay to (charm+) c+u e or fi+v^, particularly so when 
the neutrino is produced at higher energy. In r+ t~ pro- 
duction and decay, the branching ratio to v T is unity, 
but there are also 17% branching ratios to v e and to v^. 
Accompanying the is a muon, which will itself decay 
to produce another and v e , but at much lower energy 
since the muon loses considerable energy in the Sun's 



magnetic field before decaying. So roughly, the flavor ra- 
tios from the r + t~ mode arc (1.0 : 0.17 : 0.17)/1.34. In 
W + W~ production and decays, neutrinos at high ener- 
gies are democratically produced, and the associated taus 
and muons quickly lose energy in the Sun's center. How- 
ever, at lower energy, the neutrino flux originates from 
the decay of nearly-stopped heavy quarks produced in W 
fragmentation, so there is a low-energy depletion of v T s. 

A useful summary of Fig. 6 is that, at production in 
the Sun, (i) the flavor fluxes are the same for neutrino 
and antineutrino, and (ii) the e and fj, flavors are pro- 
duced in equal amounts, but the amount of r flavor de- 
pends on the WIMP annihilation channel. In the bb 
channel, r is suppressed relative to e and fi, especially 
at higher energies. In the r + t~ channel, it is just the 
opposite, with the r-flavor contribution enhanced, espe- 
cially at higher energies. The implication for the real- 
istic channel (II) consisting of 95% 66 + 5% t + t~ is a 
suppressed r component at lower energies, and a highly 
enhanced r component at higher energies. For the real- 
istic W + W~ channel (I), the r component is negligible 
below about 20 GcV, whereas t, e, and /i are equally 
represented above 20 GeV. It was claimed in Sec. II C 2 
that a democratic 1:1:1 flavor ratio at production nec- 
essarily implies a democratic flavor ratio at Earth. Thus, 
we expect flavor democracy at Earth for neutrino ener- 
gies above ~ 20 GeV, if the WIMP annihilation mode in 
the Sun is dominantly W + W~ . Figures to be shown will 
bear this out. For all other channels, and for the W + W~ 
channel below 20 GeV, the flavors at production are far 
from democratic, and the flavor-evolution machinery de- 
veloped in this work is important. 

The solar flavor spectra of Fig. 6 get processed by 
the physics we have derived in this paper. The result- 
ing flavor spectra at Earth are shown in Figs. 7-10 for 
the 6 b, t+ t~ , 95% 6 6 + 5% r+ r~ , and W+ W~ WIMP- 
annihilation channels, respectively. For each of the three 
sample WIMP masses moM = 10, 30, and 100 GeV, 
the three values 6*13 = 0°, 1°, and 12° arc considered. 
We have further taken 5 = 0, which serves to maximize 
I sin 6*i3 cos <5 1 and therefore the sensitivity to #13. Re- 
sults are displayed for the normal and the inverted mass 
hierarchies. 

Thc_ greater likelihood of the (I) W+ W~ and (II) 
95% 66 + 5% t + t~ decay channels directs our attention 
to Figs. 9 and 10. In Fig. 10 one sees the 1:1:1 flavor 
ratio at Earth for energies above 20 GeV, that results 
from the democratic production ratio in the Sun at these 
energies for the W + W~ channel. Below 20 GeV in the 
W + W~ channel, and for the whole energy range in the 
95% 66 + 5% t + channel, a comparison of Figs. 9 and 
10 to Fig. 6 makes the point that significant flavor pro- 
cessing occurs as neutrinos propagate outward through 
the Sun. The flavor processing tends to homogenize the 
\i and t components. Nevertheless, in the lower energy 
region of (I) and higher energy region of (II), the flavor 
ratios differ significantly from unity. 
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FIG. 6: Relative (anti)neutrino flavor fluxes at production in the solar core are plotted versus neutrino energy E for four 
scenarios: dark-matter annihilation predominantly into bb, or predominantly into W + W~ , or predominantly into r + r~, or 
predominantly into 95% bb + 5% t + t~. Solid (dashed) lines display the e and [i (r) flavor fluxes. For each scenario, dark- 
matter masses jtidm of 10 GeV, 30 GeV, and 100 GeV have been considered. The dark-matter mass can be inferred from the 
threshold energies of the graphs. This plot was generated from the results in Ref. [18] . 



Comparison of Figs. 9 and 10 makes the further point 
that flavor ratios arriving at Earth do vary with the 
WIMP decay channel. In fact, for the favored channels 
(I) and (II), flavor signatures are usefully in opposition, 
in that the W + W~ channel offers democratic flavors at 
higher energies and not at lower energies, whereas the 
95% 66 + 5% t + t~ channel offers nearly democratic fla- 
vors at lower energies and not at higher energies. For 
the W + W~ channel, the change from non-democratic 
to democratic is fixed at ~ 20 GeV, a number traceable 
to Standard Model masses. For the 95% 66+5% t+ t~ 
channel, the energy for change from nearly democratic to 
non-democratic depends on the WIMP mass. For either 
channel, flavor ratios in the non-democratic regions are 
seen to depend on 6*13, and, for 6*13 ^ 0, on the mass hier- 
archy. For the 95% 66 + 5% t + t~ channel, certain flavor 
ratios may be as large as 2 at the higher energies; and 
in the W + W~ channel, as large as 1.4 below ~ 20 GeV. 
The flavor ratio of 2 is within experimental reach, while 
ratios of 1.4 or smaller may present a challenge for some 
neutrino telescopes [21] but perhaps not for others [22]. 

IV. SUMMARY AND CONCLUSIONS 

The discovery of > GeV neutrinos from the Sun would 
constitute strong evidence for the annihilation of solar 



WIMPs. Competing WIMP models can in principle be 
distinguished by the energy spectrum and flavor ratios 
of the decay neutrinos and antineutrinos. The flavor ra- 
tios detected at Earth will differ from the flavor ratios 
at injection in the Sun's core. Due to the decoherence 
resulting from ( L/_E)-averaging in the Sun, the flavor 
evolution is described by classical probabilities, rather 
than by quantum-mechanical amplitudes. This simplifies 
the analysis. On the other hand, the evolution is com- 
plicated by the matter effects within the Sun. WIMPS 
annihilate in a CP-symmetric state, and so the neutrino 
and antincutrino spectra at the point of origin are iden- 
tical in energies and flavors. However, since the material 
of the Sun is matter and not antimatter, propagation 
effects through the Sun are different for neutrino ver- 
sus antincutrino. In the travel outward from the Sun's 
core, two level-crossings are encountered (one if #13 = 0), 
each described by an MSW resonance. For the normal 
neutrino-mass hierarchy, both level crossings occur in the 
neutrino sector; whereas for the inverted hierarchy, the 
lower-energy crossing occurs in the neutrino sector and 
the higher-energy resonance occurs in the antineutrino 
sector. Thus, the emergent energy and flavor spectra de- 
pend also on the neutrino-mass hierarchy, as indicated 
schematically in Fig. 1. An additional complication is 
that evolution at the level crossings may be adiabatic or 
non-adiabatic, depending on solar and neutrino parame- 
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FIG. 7: Terrestrial fluxes versus neutrino energy E for solar WIMP annihilation predominantly into bb. The quantity (23V + 
$-)/3 is plotted, where the v and v fluxes are weighted according to their cross sections for detection. This quantity would 
be measured by Earth-based detectors that cannot resolve v a from V a events, as discussed in Sec. II B 5. The neutrino mass- 
squared differences are <5m!i = 8.0 x 10 -5 eV 2 and 8m%2 = ±3.0 x 10~ 3 eV 2 . The mixing angles are 812 = 35° and 623 = 45°, 
and the CP- violating phase 5 has been set to zero. We have considered the values 813 — 0° , 1°, and 12° for each hierarchy and 
for each of the dark- matter masses hidm = lOGeV, 30GeV, and 100 GeV (from top to bottom). The and v T spectra as well 
as the Ufi and u T spectra are interchanged if 8 — ir is chosen. Note that the observables depend noticeably on the mixing angle 
#13. Moreover, for larger 813 angles there are visible differences between the normal and inverted hierarchies. 



ters. 

In this work, we presented analytic formulas for the 
complete flavor evolution. One virtue of an analytic 
treatment is the emergence of post hoc intuition. With 
the gift of hindsight, we therefore were able to present 
the physics of flavor evolution which underlies the 
analytic formulas. For example, we were able to derive 
general statements that 

(i) If flavor democracy, defined by 1 : 1 : 1 flavor ratios, 
occurs in WIMP decay, then the same 1:1:1 ratio must 
be observed on Earth; in other words, flavor-mixing, 
even with complicated matter effects, cannot derail 
democracy! 

(ii) For maximal 6*32 but arbitrary #21 and 9±3, the 
replacement 5 —* ir — S interchanges the role of and 
v T , and and P r , and leaves unaffected the terrestrial 
v e and v e fluxes. 

(iii) For #32 and 6*12 at their tribimaximal values, the 
terrestrial flavor ratios for normal-hierarchy neutrinos 
and inverted-hierarchy antineutrinos are independent 
of 5 to zeroth order in (Sm^i/Sm^)- In addition, the 
terrestrial v e and v e fluxes are independent of 5 to 
zeroth order in (Sm^/Sm^) for any hierarchy. (This 
feature holds only in a finite energy region determined 
by adiabaticity and by V e E ^> Sm^ k .) 

Concerning the adiabaticity versus non-adiabaticity, 
we found that level-crossings are adiabatic in the neu- 



trino sector for energies up to ~ 10 GeV or more in the 
neutrino sector. In the antineutrino sector, there are no 
level crossing for the normal hierarchy and so adiabatic- 
ity holds trivially at all energies, while for the inverted 
hierarchy, adiabaticity holds up to <~ 50 GeV or more. 
The energy-onsets of non-adiabatic behavior (En a) are 
indicated in Figs. 2-4. 

There is a striking dependence on the value of the un- 
known mixing-angle #13 in the evolution of flavor away 
from a single injected flavor in the Sun. In particular, the 
limit 9 1 3 — > is "discontinuous" in that for 13 = the 
third neutrino mass eigenstate | 3 ) is completely decou- 
pled from the other two mass eigenstates, and there is but 
a single resonance is encountered in the flavor evolution; 
while for 6*13 7^ 0, two resonances are encountered, and 
even a small nonzero value of #13 may be amplified by 
the associated resonance. In addition, for nonzero #13, 
there is sensitivity to the CP- and T-symmctric term 
sin 6*i3 cos S, where S is the CP- and T- violating phase 
parameter in the unitary mixing-matrix. Flavor evolu- 
tion from a single injected flavor in the Sun, for #13 set 
to zero, is shown in Fig. 2. For the maximal allowed 
value of 9is , the same is shown in Fig. 3 for the normal 
mass-hierarchy, and in Fig. 4 for the inverted hierarchy. 

We further explored the dependence on S in Fig. 5. 
The dependence on 6 enters through the CP- and T-even 
factor 3?{[/ e 3} = sin #13 cos<5. To maximize sensitivity to 
6 in this figure, we have again chosen for 613 the max- 
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FIG. 8: Terrestrial fluxes versus neutrino energy E for solar WIMP annihilation predominantly into t + t~. Again, the observable 
(2$,, + )/3 is plotted. Neutrino parameters are identical to those employed in Fig. 7. The and v T spectra as well as the 
Vfi and v T spectra are interchanged if S — tt is chosen. Sensitivities to 813 and the mass hierarchy occur at levels similar to 
those of the bb case in Fig. 7. 
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FIG. 9: Terrestrial fluxes versus neutrino energy E for solar WIMP annihilation into 95% bb + 5% t + t~. Again the observable 
(23v + &v)/3 is plotted, and neutrino parameters are identical to those in Figs. 7 and 8. The and v T spectra as well as the 
and v T spectra are interchanged if 8 = tv is chosen. Sensitivities to #13 and the mass hierarchy are comparable to those in 
Figs. 7 and 8, as expected. 



imum value allowed by the CHOOZ data, 12°. The jj, 
and t flavored spectra are sensitive to 5, whereas the e 
spectrum is not. The and v T spectra and the and 
v T spectra swap among themselves for 5 — > n — S. These 
results are in accord with the general features mentioned 



above. 

WIMP decay does not produce neutrinos of just one 
flavor. Consequently, our illustrative evolution results in 
Figs. 2-5 do not apply directly to the flavor processing 
of WIMP-decay neutrinos. Rather, some linear com- 
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FIG. 10: Terrestrial fluxes versus neutrino energy E for WIMP annihilation predominantly into W + W~ . The plotted observable 
(2<I?„ + <&v)/3 and the neutrino parameters are the same as those in Figs. 7-9. The and v T spectra as well as the and v T 
spectra are interchanged if S — ir is chosen. Sensitivities to #13 and the mass hierarchy are comparable to those in Figs. 7-9. 



binations of the single- flavor figures are relevant, fn 
most WIMP models, the annihilation final states fall 
into one of two categories. The two are the W + W~ fi- 
nal state and the 95%, 66 + 5%t+t~ final state. These 
Standard Model final states decay further in well known 
ways to produce the neutrinos and antineutrinos ema- 
nating from the Sun's core. In Fig. 6 we showed the 
injection neutrino and antineutrino spectra per flavor, 
gleaned from [18]. A representative sample of WIMP 
masses were chosen: 10 GeV, 30 GeV, and 100 GeV. 
The flavor spectra arriving at Earth, after mixing and 
matter processing, we showed in Figs. 7-10 for various 
possible WIMP-annihilation final states. Each of these 
Figures is arranged into panels labeled by the three val- 
ues 0i3 = 0°, 1°, 12°, and the three values of WIMP mass 
10 GeV, 30 GeV, 100 GeV. The phase 5 has been set to 
zero (or 5 = it with the interchange of and v T spec- 
tra and and v T spectra); we reserve for future work 
the dependence of the observable flavor spectra on 6 [35]. 
One feature of the flavor evolution is a tendency to av- 
erage the \i and r flavors. Nevertheless, some interesting 
nontrivial flavor ratios appear in these figures. 

For the two preferred channels of WIMP annihila- 
tion, it is seen that in the lower energy region of the 
W + W~ decay mode, and in the higher energy region of 
the 95% 66+5% r + t~ decay mode, flavor ratios differ 
significantly from unity. Conversely, in the higher and 
lower energy regions, respectively, the flavor ratios are 
(nearly) democratic. Since the decay mode implicates 
the nature of the WIMP, the experimental ability to dis- 
tinguish decay modes allows an inference of the nature of 
the WIMP. The energy dependences of the neutrino fla- 
vor ratios at Earth arc, then, the signatures which poten- 
tially distinguish the two main competing WIMP models. 

The flavor ratios at Earth depend on #13 and on the 
neutrino mass hierarchy. For the purpose of studying the 
nature of solar WIMPS, it would be very useful if near 
future experiments would better determine these two pa- 
rameters. On the other hand, if the nature of the WIMP 
becomes known from future accelerators such as the LHC 
and/or the ILC, then it becomes conceivable to use the 
flavor ratios to deduce or bound better 0i 3 , and possibly 
deduce the mass hierarchy. 



If one attempts either to infer the nature of the so- 
lar WIMP, or to determine the values of the parameters 
013 and sign(<5m| 2 ), an algorithm for neutrino flavor pro- 
cessing from the Sun's core outward is essential. In this 
paper, we have provided this algorithm. 
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APPENDIX A: ANALYTIC U m FOR NONZERO 
0i3 AND S 

In order to obtain a manageable expression for U m , we 
restrict our attention primarily to perturbations about 
the tribimaximal case (10). In particular, we take 
012 = 7r/6 + (50i2, and we implicitly assume tan0i3 <, 
O(l), which is observationally justified. On phenomeno- 
logical grounds, we also have <5m|i/ '8m\ 2 "C 1 and 
Srrij k /V e (0)E < 1, so that these quantities can serve as 
additional expansion parameters. In the latter relation, 
we have employed the fact that our present focus is on 
GeV neutrinos. 

We begin by recalling the defining relation for U m : 



H M = Ul UMU^ + V U, 



(Al) 



such that Hm is diagonal. For notational brevity we have 
defined M = diag(— (Jm^i, 0, 5m\^)/2E (as before), and 
V = V e diag(l, 0, 0). In order to simplify the algebra, we 
decompose U m as follows: 



u m = R 23 (e 2 3)ulu n 



(A2) 
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with U m to be determined. Employing this decomposi- 
tion and Eq. (7) in the definition (Al) we obtain 



H M (5) = ill 



R 13 U S R 12 M R\ 2 ujR\ 3 + V 



Um - (A3) 



Here, we have suppressed the arguments of the rota- 
tion matrices, and we have used that both R23 and Us 
commute with V. Since U m is unitary and diagonal- 
izes Ri 3 UsRi2 M R\ 2 UlR.\ 3 + V = H, it is composed 
of normalized eigenvectors of H. To find these, we can 
employ perturbation theory: we have 5m 2 k <C EV e , so 

that R 13 U s Ri2 M R\ 2 uIr\ 3 = SH can be treated as a 
perturbation when compared to V. 



The first step in time-independent perturbation theory 
is to identify the zeroth-order eigenstates, which in the 
present situation are those of V. Up to a phase, the first 
such eigenvector clearly is 



vi = (1,0,0) . 



(A4) 



The remaining two eigenvectors v 2 and v 3 , which are a 
linear combination of (0, 1, 0) and (0, 0, 1), are associated 
with degenerate eigenvalues. They must therefore be de- 
termined by diagonalizing SH restricted to the 2-3 sub- 
space. This essentially amounts to finding the eigenvec- 
tors of the lower right 2x2 block of 



SH 



1 



6 Sm 32 sin 2 i3 — 4 Sm 21 cos 2 9\ 3 
— , 2V2e-* 5 / 2 5m 2 21 cos 9 13 
(2 Sm\ x + 3 Sml 2 ) sin20i 3 

° 2a ~- e u / 2 cos9 13 



2V2 e lS/2 Sm 2 21 cos 6 13 (2 Sm\ x + 3 Smj 2 ) sin 26 13 
-2Sm 2 21 -2^2 er l5 / 2 Sm 2 21 sin 13 

-2^/2 e iS / 2 Sm^ sm6 13 6 Sm\ 2 cos 2 9 13 - 4 Sm\ x sin 2 6 13 

-\/2 sin 26>i a \ 



3Sm 32 )sin29 13 



c 2 /2V2cos 2 13 e l °^a 
+ e-^/ 2 cose 13 -2V2 



-s/2 sin 26»i3 

, . — ■ \ — -e~ iS / 2 sin 6<i3 | <S0i 2 

\-V2sin26»i3 -e l5 / 2 sin 6»i3 2V2sin 2 0i 3 



(A5) 



Here, 8O12 denotes deviations of Q12 from its tribimaximal 
value implicitly determined by sin 2 812 = 1/3. The above 
expression (A5) is valid at first order in 6612. We remark 
that the correction term involving 58i 2 contains Sm 21 as 
an overall factor, but Sm 32 is absent in the SQ\ 2 term. 
This essentially arises because in the definition of SH the 
i?i2 matrix mixes only the upper 2x2 block of M, so 
that any occurrence of #12 is necessarily accompanied by 
a Sm 21 factor. 

For the remaining zeroth-order eigenvectors, we now 
find at leading order in Sm 21 /Sm 32 



( 



V2 





e -iS/2 

\/2 tan 013 Sm 21 



\ 



(A6) 



\ 3 cos6>i3 5m\ 2 / 



and 



v 3 = 



\/2 tan 0i3 Sm 21 
3 cos 0i3 Sm 2 



iS/2 



32 



V 



I 



(A7) 



We remind the reader that normalized eigenvectors are 
only defined up to phases. Note that these expressions 
arc independent of <50i2- This is a consequence of the 
fact that Sm 32 is absent in the S9i 2 term in Eq. (A5), so 
that the <50i2 correction contains a further suppression 
relative to the <50i2-independcnt contribution in SH. 



The corrections to these eigenvectors are given by the 
usual perturbation-theory formula 



Svj 



k SHv J 
Xj — Xk 



Vk 



(A8) 



where Xj denotes the eigenvalue of V corresponding to 
Vj, i.e., Ai = V e and A 2 = A3 = 0. The sum in Eq. 
(A8) runs over all eigenvectors that do not belong to the 
eigenspace containing vj. We obtain 



/ 



Svi 



Sv 2 = 



5v 3 = 



V e E 



V e E 



V R E 







^e^ 5 ' 2 Sm 2 21 cos 9 13 



V ]^(3 5m| 2 + 2<5m^)sin20i 3 / 
1 

3V2 



( — 5m 21 cos 0i3 \ 



V 



/ - ± (3 Sm\ 2 + 2 5m 2 2l ) sin 20 i3 \ 





(A9) 



at leading order in Sm 21 /5m\ 2 . It follows that our ap- 
proximation U m ~ (vl + Svi , v2 + Sv 2 , v3 + Sv 3 ) is explic- 
itly given by 
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U„ 



Sm 21 e 



-iS/2 



COS ^13 



3V2V e E 
3 Jm 2 ^) sin 20 



(2Sm 2 21 



13 



12V e £ 



<5m 2 icos0i 3 (2 (5m2! + 3(5m§ 2 ) sin2^i3 \ 



3^1^ 

e -M/2 

^/2Sm 21 tan 

3(5TO3 2 COS 013 



12V e E 
V2e- lS ^ 2 8m 2 21 tan gig 

3 #77232 COS 6*13 
1 



(A10) 



which is valid to linear order in both 8m 2 2 /V e E and 
8m2i/5m\ 2 . First-order terms in <50i 2 are absent. This 
again arises because the £0i 2 contributions in Eq. (A5) 
suffer from an additional suppression relative to the 
leading-order terms in 5H. As pointed out in Sec. II B, 
the results for antineutrinos can be obtained by reversing 
the signs of V e and 5 in the above equations. 

Note that in U m and thus in U m , corrections contain- 
ing 0i3 or S are also suppressed by Sm 21 /5m 2 2 <C 1 or 
5m 2 k /V e (0)E <C 1. This implies that for GeV neutrinos 
at the solar core, the matter mixing is largely unaffected 
by 0i3 and d perturbations about the tribimaximal case. 
Significant 0i 3 and S effects on the terrestrial fluxes can 



therefore only arise through the non-adiabatic transitions 
described by the matrix P and through the vacuum mix- 
ing matrix U, which determines the flavor content of the 
mass states at the detector. 

We remind the reader that each column vector in U m 
is only defined up to a phase, an ambiguity without phys- 
ical significance. The column- vector order in U m , on the 
other hand, is to be selected such that the diagonal el- 
ements of U^HUm exhibit the same ordering as those 
of the matrix M. The arbitrary ordering choice in Eq. 
(A10) gives 



UlHU m = — 



2V e E + 5ml 2 sin 2 13 







i + • • ■ 

Sm 2 2 cos 2 0i3 



(All) 



where the ellipsis indicates terms of second or higher or- 
der in Sm 2 2 /V e E and 5m 21 /Sm 2 2 . The present range of 
parameters yields <C |3cos 2 0i35m 2 2 \ <C |6V e _E|, where 
the first inequality rests on our assumption that 0i 3 is 
not too large. The correct ordering therefore depends on 
the signs of V e and 8m\ 2 . It follows that the required 
adjustments can only be implemented once particle type 



(y vs. v) and mass hierarchy are specified. 

Apart from the column- vector order in U m , Eq. (A2) 
can be used to find explicit first-order expressions for U m , 
and thus U m . The matrix U m takes the relatively simple 
form 



U„ 



1 

cos 2 23 sin 2 23 | - 
Sin 2 023 cos 2 23 



cos S sin 0i3 sin 20 2 3 

3 COS 2 013 




1 

-1 



5ml 2 



(A12) 



This expression does not contain 0(Sm 2 k /V e E) and 
0(6612) terms. Only the sign of V e will eventually play 
a role, as it matters for the column-vector ordering in 
U m . It is again apparent that the effects of 0i3 and 6 on 
the matrix U m are suppressed by more than an order of 
magnitude from the 5m 2l /5m\ 2 ratio. For the observ- 
able T-conserving matrix elements of U m , the 6 depen- 



dence can only enter in the combination ^R.(e tS sin0i 3 ) = 
cos 6 sin 0i 3 . 
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APPENDIX B: EQUALITY OF LOW- AND 
HIGH-ENERGY FLAVOR RATIOS FOR i3 = 







In the low-energy limit E — > 0, matter effects are neg- 
ligible, so that PU^ — > ?7 T in Eq. (44). In other words, 
the vacuum case (18) arises. This general result must 
hold for all values of the #13 angle. In the high-energy 
limit E — > 00, both P and U m in Eq. (44) are nontrivial 
because all matter effects must be taken into account. It 
follows that the low- and high-energy limits are distin- 
guished by the difference in the U T and PUj n matrices, 
where the latter is to be evaluated for neutrinos (as op- 
posed to antineutrinos) in the limit E — > 00. We will 
show that U T and PU^ are identical when #13 = 0. 

For #13 — 0, the Hamiltonian can be transformed into 
the form H' F given by Eq. (24). As noted in the discus- 
sion leading to Eq. (29), the off-diagonal pieces of H' F 
vanish when E is taken to approach infinity. It follows 
that U m = i?23(#23) is one possible matrix that diagonal- 
izes the original Hamiltonian Hp at high energies. The 
remaining ambiguity (up to phases) in the construction of 
U m is determined with our column-vector ordering con- 
vention. To this end, we need to fix the mass hierarchy 
and take 5m\ 2 > as an example. The reader is invited 
to verify that the inverted hierarchy produces the same 
final result, as expected due to the decoupling of the |3, r) 
state for 6*13 = 0. Note that V^O) is positive because are 
considering neutrinos, and not antineutrinos. It follows 
that the largest eigenvalue in H' F is given by V^O) lo- 
cated in the upper left entry. To match our convention, 
this eigenvalue should appear as bottom right matrix el- 
ement; the remaining two eigenvalues are already in the 
correct order. We thus modify our U m above with an 
appropriate permutation matrix, to get 



U m = R 23 (8 23 ) 




(Bl) 



as the transformation that diagonalizes the high-energy 
flavor-basis Hamiltonian for 6*13 = 0. 

Next, we turn to the crossing-probability matrix P. 
For E > E NA , P l c -> cos 2 6>i2 and P c h -» cos 2 ff 13 = 1 
at the respective lower and higher resonances. Note that 
the latter of these probabilities is consistent with the fact 
that the |3,r) state decouples for #13 = 0. Employing 
these results in Eq. (55) leads to 



sin 2 612 cos 2 
cos 2 612 sin 2 9 l2 





(B2) 



for the crossing probabilities at high energies in the case 
of a vanishing #13 angle. 

We are now in the position to give the simple expres- 
sion 



/cos 2 #i2 sin 2 #i2 cos 2 #23 sin 2 #i2 sin 2 # 2 3 



PU T 



sm 2 #i2 cos 2 #i2 cos 2 #23 cos 2 $i 2 sin 2 f 



23 



sin 2 (9 



23 



cos 2 #- 
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for the high-energy limit of the matrix PU^ for neutri- 
nos. Comparison of this result with U T then establishes 
that the these two matrices are identical. Therefore, 
the low- and high-energy terrestrial neutrino flavor ratios 
are also identical. Moreover, the low-energy antineutrino 
fluxes, which are obtained by replacing U — > U* in Eq. 
(17), are equal to the above low-energy neutrino results. 
Both of these features are apparent in Fig. 2. 



APPENDIX C: PROOF OF FEATURE (i) OF 
SEC. IIC 6 

We first note that the level-crossing probabilities pre- 
sented in Sec. II C 2 act on mass eigenstates, and so are 
independent of S. It follows that <5 can only affect the ter- 
restrial neutrino fluxes via the U and U m flavor-content 
matrices. This implies that possible phase effects may 
enter only at the production and detection sites. That 
neutrino propagation including non-adiabatic transitions 
in the Sun are unaffected by 5 does not come as a surprise, 
as our jump-probability approximation involves only two 
levels at each resonance, and in two-flavor systems the 
phase in the mixing matrix can be removed. 

Thus, to show that the replacement 6 — > ir — S leaves 
unaffected the v e and v e spectra, while interchanging the 
Up and v T spectra, and and v T spectra, it is sufficient 
to show that the replacement leaves unaffected the first 
row in both U and U m , but interchanges the lower two 
rows. In mathematical terms, we must show that 



(U)ej(S) = (JDej^-S), 

(U)M = (U) TJ (n-5), 



and 



{Um)eM 



(U m ) Tj (ir - 5) , 



for.? = 1,2, 3. 

Explicit evaluation of U (S) and U (tt—6) with # 32 
and 621 and #13 arbitrary yields 



(CI) 



(C2) 



45° 



U ej (5) 
U ej (6) 



U ej (n-S) 
-U* ej (-K-S) 
-U* Tj {-K-5) 
U Tj (n-S) 



J = 1,2 
3 = 3 , 
J = 1,2 
3 = 3 • 



(C3) 
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These relations imply our condition (CI), as can be ver- 
ified by multiplying each equation with its complex con- 
jugate. We remark that this result already establishes 
Feature (i) for cases with negligible matter effects, such 
as in vacuum. 

Large matter potentials V e (r) can definitely modify 
neutrino flavor fluxes, as we have seen in the text. How- 
ever, the Vfj,-v T relationship is left unaffected by V e (r) 
because only v e interacts with the solar background elec- 
trons. The Vn-v T part of the above vacuum result, 



27 



namely that the and v T fluxes switch under 5 — > it — 5, 
should therefore generalize to situations with matter. 

To see this explicitly and to establish furthermore that 
the v e flux remains unchanged, we need to prove condi- 
tion (C2). The mixing matrix in matter is defined via Eq. 
(Al). With the 5 dependence displayed, this definition 
reads for U m (S): 

H M (S) = U m (Sy [U(6)MU\5) + V] U m (S) . (C4) 

In order to relate U m (n — S) and U m (5), we start from 
the analogous equation for U m (ir — S), i.e., Hm(k — 8) = 

W m (it - S) \lf(ir - S) M W(n - S) + v) U m (ir - S), and 

implement the symmetries (C3). In matrix form, these 
symmetries may be rewritten as 



U(ir-S) = FU*(6)S , 



(C5) 



where the matrix F interchanges the /i and r flavors, and 
the matrices S and F change various signs: 



S = 



-1 
















1 


) 





1 







-1 















1 













1 



(C6) 



The specific expressions for the matrices M, V, F, and 
S are such that SMS^ = M and FVF^ = V. With 
these observations, the defining relation for U m (ir — 8) 
can be cast into the following form: Hm(v: — 8) = 

WJt:-5)f\u*(8)MU t (5) + V 



F^U m (w-S). Recall- 
= F* are real, complex 



ing that H M (n-5), M, V, and F = 
conjugation of this equation yields: 



Hm(tt-S) = 

E£(tt -5)F [U(6) M UHS) + V] Ft U* m {-K - 8) . (C7) 

By definition, both Hm (7r — 8) and Hm (8) are diagonal. 
Moreover, they exhibit the same ordering of eigenvalues 
because they are associated with the same mass hierarchy 
and particle type [y vs. 77). Note also that F^ U^tt — S) 
is unitary. These facts and comparison of Eq. (C7) with 
Eq. (C4) show that F^ U* n (ir - 5) satisfies the defining 
relation (C4) for U m (8). We thus find 

/ e M 

U m (6)= | F^U* m {ir-5). (C8) 

\ e v 

Here, the three undetermined phases a, /?, and 7 arc 
unobservablc, but are included for completeness. The 
condition (C2) follows now directly from (C8). Feature 
(i) is thereby established. 



APPENDIX D: PROOF OF FEATURE (ii) OF 
SEC. IIC 6 



The general terrestrial flavor fluxes are given by Eq. 
(44). The matrix A = UPUj n in this equation can 
explicitly be obtained at first order with the result 
(A12). At leading order in the three small parameters 
Sm^x/Sm^j S9i2, and <5#23, the generic form of A is 



A = A + A, 



8n4i 
'■ 5mj 2 



+ A 12 (56»i2 + A 23 69 23 



(Dl) 



Here, S9i 2 and S623 denote deviations from the tribimax- 
imal values of 6*12 and 623, respectively. Note that terms 
of order 5m^ k /V e E are absent because they do not enter 
into U m in Eq. (A12), and U in vacuum depends only 
on the mixing angles. Below we give explicit expressions 
for the matrices Ao , A m , A12, and ^23 • These explicit 
expressions require the knowledge of the column-vector 
ordering in Eq. (A12). Thus, the specification of particle 
type (y vs. V) and the mass hierarchy is necessary. 

We seek to prove that in the adiabatic limit (P = 1), 
for NH neutrinos and IH antineutrinos, that Ao is inde- 
pendent of S; and that with either mass hierarchy, the 
first row of A is independent of S. When A a is inde- 
pendent of (5, then all flavor spectra are independent of 
S at leading order; when the first row of A n is indepen- 
dent of S, then the v e and v e spectra are independent of 
S at leading order. Below we determine the full 6*13 and 6 
dependence of A . For completeness, we also determine 
A12 and A23 at zeroth order in the small parameters. It 
turns out that A m = + 0(813, 59\ 2l S823) in all cases, 
and so we do not consider it any further. 

To estimate the effects of non-adiabaticity (P 7^ 1), 
we also consider the case of a completely non-adiabatic 
higher-energy resonance while maintaining adiabaticity 
at the lower-energy resonance. This requires that #13 
be small enough such that the non-adiabaticity at the 
higher-energy resonance occurs at energies where the 
lower-energy resonance is still adiabatic. 

Neutrinos within the normal hierarchy. This case 
is characterized by < 8m\ 2 < V e E. Inspection of Eq. 
(All) shows that the U m column vectors in Eq. (A10) 
need to be rearranged according tol— »3— »2— >1, i.e., 




(D2) 



This gives the following explicit form of the A , A 12 , 
A23 matrices: 



and 



An = 



4s 2 

zc 13 
Z6 13 



zc 13 



2r 2 

ZL 13 



1 + ,S?. 



1 + s 



'13 
2 
13 



A 12 = , 



(D3) 



(D4) 
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-423 = TT 




(D5) 



Inspection of Eq. (D3) shows that S leaves unaffected the 
terrestrial flavor fluxes at leading order, which establishes 
Feature (ii) for normal-hierarchy neutrinos. Moreover, 
the second and third rows of A are identical, implying 
degeneracy between the and v T fluxes at this order. 
The results (D3), (D4), and (D5) become invalid for 



small #13. In this situation, the resonance involving |3, r) 
fails to meet the above adiabaticity condition P = 1: the 
onset of non-adiabaticity (49) applies here with 6 = 6*13, 
which lies below the resonance energy (53) for 6*13 small 
enough. This non-adiabaticity is associated with the de- 
coupling of the state |3,r), as discussed in Sec. IIB1. We 
display the A-matrices for a completely adiabatic lower- 
energy resonance and a completely non-adiabatic higher- 
energy resonance: 



/c 13 



2 + S13 ~ \/8C5Sl3 2~\s\ 3 
\2 + s 2 13 + V8css 13 2-\s\ 3 



2 + s 2 3 



\ 



-\/2c«5Si3 2 
V2 ess 13 2 



is 2 
2*13 

2 S 13 



- V%CS S13 
V2csSi 3 ) 



(D6) 




In this non-adiabatic case, only the first row of Aq is 
independent of 5, and so only the v e and D e spectra are 



independent of 5 at lowest order. 

Antineutrinos within the normal hierarchy. As 

the sign of V e for antineutrinos is opposite that for neu- 
trinos, we have now V e E < < <5m§ 2 . It follows that the 
diagonal entries of U m HU m in Eq. (All) are correctly 
ordered already, and adjustments in U m are unnecessary. 
But for antineutrinos we need and U* , which can be 
obtained by simply changing the sign of S in the final 
expressions for A , A12, and A 23 . This yields 



A = 



4c? 



13 1 + 2s? 3 1 + 2s 2 3 

V8c s s 13 ^ + cj 3 -V2css 13 § + c\ 3 - y/2cssi 3 
\-Z-2c\ 3 - V8c s s 13 %+ci 3 + V2c s s 13 | + c 2 3 + V2c s s 13 J 



3 - 2c\ 3 



(D9) 




(DIO) 



1 ( -1 1\ 
A 23 = - -1 1 . (Dll) 
6 \ 1 -1 J 

Note that there arc unsupprcssed phase effects, since Aq 
depends on S. However, the first row of Aq remains 6- 
independent. 

We remind the reader that resonances are absent in 
this case. Adiabaticity is therefore ensured trivially, and 
the results (D9), (DIO), and (Dll) remain valid in the 
limit 613 -> 0. 



Neutrinos within the inverted hierarchy. This 
case is characterized by 5m 32 < < V e E. Inspection 
of Eq. (All) and M establishes that we have to inter- 
change the first and second column in U m , which can be 
implemented via the replacement 



Um -> U m 1 . (D12) 



We then obtain the following expressions for A , Ai 2 , 
and A 23 : 
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An = 



2 — 2sf 3 
2 + S13 - V&c s s 13 



2 + s 2 



13 



2 + s 2 



13 



\ 



2- 



|s? 3 + V2c 5 si3 2 - |sf 3 + \/2c 5 si3 
V 2 + s 2 3 + y/8c s S13 2 - isf 3 - V2 c 5 S13 2 - ±sf 3 - V2c 5 S13 y 



(D13) 



! / 4 -2 -2\ 

-**t:5 ; ;)' (m4) 



for neutrinos within the normal hierarchy: 



= 7 



4 2 or 2 2r 2 
^*13 zc 13 z6 13 
„2 1 , „2 1 , „2 



2cf 3 1 + S f 3 1 + S f 3 | , (D17) 

2 C 13 1 + s 13 1 + s 13 




(D15) 



It is apparent that phase effects are again unsuppressed 
because An depends on S. However, the first row of An 
remains <5-independent. 

Although a resonance is present in the neutrino sector 
it does not involve the state |3, r). Hence, the adiabatic- 
ity condition is met, and the results (D13), (D14), and 
(D15) keep their validity for vanishing 9\ 3 . 

Antineutrinos within the inverted hierarchy. In 
this case, we have V e E < <5m 32 < 0. Comparison 
of the eigenvalue ordering in Eq. (All) with that in 
M = diag(— 5m 2 ;!, 0, Sm 32 )/^E shows that we have to 
interchange the first and third column in U m , which can 
be implemented via the replacement 

- ~ (° ° M 

U m - U m 1 . (D16) 

\1 0) 

The final results for An, A\ 2 , and A23 are similar to those 



A 12 = , (D18) 



1 ( ° 1 ~ X \ 

A 23 = - 3 -2 -1 . (D19) 

6 \ -3 1 2 J 

We see that the present An is identical to that in Eq. 
(D3). It follows that within the adiabatic regime, the 
net behavior of IH antineutrinos equals that of NH neu- 
trinos within the approximations made in this appendix. 
Feature (ii) therefore also holds in the present case. 

Let us mention an additional aspect shared by NH neu- 
trinos and IH antineutrinos: the presence of the reso- 
nance involving |3,r). Hence, the results (D17), (D18), 
and (D19) become invalid for small 6*i 3 , paralleling the 
above argument for NH neutrinos. Finally, we quote 
the result for the completely non-adiabatic limit at the 
higher-energy resonance: 



An = - 



/ 

3 

\3 



4c 2 



2c\ 3 

2r 2 
zc 13 



13 

VScg S13 



1 + 2s? 



13 

§ + c 2 3 - V2c s s 13 



1 + r 2 

2 T C 



1 + 2s 2 3 

V2C5 S13 



13 



(D20) 



V8c s s 13 § + c 2 13 + V2c s s 13 | + c 2 3 + V2c s s 13 J 




Summary — To summarize this appendix, we have 
shown that when 9 32 and #21 assume their tribimax- 
imal values, then to zeroth order in Sm^/Sm^ and 



6m% 2 /V e (0) E, one has that 

(i) the v e and v e spectra are independent of the phase- 
parameter 8; 

(ii) thc normal-hierarchy and v T spectra are indepen- 
dent of 5; and 

(iii) the inverted-hierarchy and v T spectra are inde- 
pendent of 5. 

Only the normal-hierarchy and v T spectra, and the 
inverted-hierarchy and v T spectra depend on S in ze- 
roth order. Finally, we remark that the dependence on 
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S in the derived expressions has the form sin #13 cos<5, as required for CP- and T-invariant observables. 
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